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Abstract 



The electroweak phase transition is investigated by means of the perturbatively calculated 
high temperature effective potential. An analytic result to order g 4 , A 2 is presented for 
the Abelian Higgs model, the SU(2)-Higgs model and the standard model and a complete 
on-shell renormalization at zero temperature is performed. Higher order corrections are 
found to increase the strength of the first order phase transition in the non- Abelian model, 
opposite to the Abelian case. This effect is traced back to the infrared contributions from 
the typical non- Abelian diagrams. The dependence of several phase transition parameters 
on the Higgs mass is analysed in detail. A new, gauge invariant, approach based on the 
composite field is introduced. This method, which supports the above Landau gauge 
results numerically, permits a conceptually simpler treatment of the thermodynamics of the 
phase transition. In particular, it enables a straightforward comparison with lattice data 
and the application of the Clausius-Clapeyron equation to the electroweak phase transition. 



Zusammenfassung 

Der elektroschwache Phasenubergang wird mit Hilfe des storungstheoretisch berechneten 
effektiven Potentials der Hochtemperaturfeldtheorie untersucht. Fur das abelsche Higgs- 
Modell, das SU(2)-Higgs-Modell und das Standardmodell werden analytische Ergebnisse der 
Ordnung g 4 , A 2 angegeben, und eine vollstandige Renormierung im On-Shell-Schema bei 
Temperatur Null wird durchgefuhrt. Im Gegensatz zum abelschen Fall fiihren die Korrektu- 
ren hoherer Ordnung beim nichtabelschen Modell zu einer Verstarkung des Phasenubergan- 
ges erster Ordnung. Dieser Effekt wird auf die infraroten Beitrage typischer nichtabelscher 
Diagramme zuruckgefuhrt. Die Abhangigkeit mehrerer Parameter des Phasenuberganges 
von der Higgs-Masse wird einer detaillierten Analyse unterzogen. Ein neuer, eichinvarianter 
Zugang, der sich auf das zusammengesetzte Feld stiitzt, wird eingefuhrt. Dicse Metho- 
de bestatigt numerisch die obigen, in Landau-Eichung erzielten Resultate und erlaubt eine 
begrifflich einfachere Behandlung der Thermodynamik des Phasenuberganges. Insbeson- 
dere ermoglicht sie den unmittelbaren Vergleich mit Gitter-Daten und die Anwendung der 
Clausius-Clapeyron- Gleichung auf den elektroschwachen Phasenubergang. 
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Introduction 



The whole known universe consists almost exclusively of matter, with no considerable amount 
of antimatter in our galaxy cluster and no known mechanism to separate matter and anti- 
matter on such large scales [ffj. This baryon asymmetry of the universe is one of the most 
interesting cosmological problems to be resolved by particle theory. 

As pointed out by Sakharov in 1967 the baryon asymmetry may be a calculable result 
of particle interactions || , the necessary conditions being baryon number violation, C- and 
CP-violation, and departure from thermal equilibrium. 

Kirzhnits and Linde realized in 1972 that at high temperatures the spontaneously broken 
electroweak symmetry is restored, thus suggesting a phase transition in the early universe 
0. Since anomalous baryon number violation in the standard model is rapid at high tem- 
peratures, the departure from equilibrium in a first order electroweak phase transition opens 
the possibility of standard model baryogenesis. This scenario, first suggested by Kuzmin, 
Rubakov and Shaposhnikov in 1985 |4]], provides the main motivation for the present investi- 
gation. Although more recent analyses seem to discourage baryogenesis within the minimal 
standard model due to the small CP-violation, simple non-minimal models may produce a 
sufficient asymmetry (see || and references therein). In any case it is clear that the present 
baryon asymmetry of the universe has been finally determined at the electroweak phase 
transition, since baryon-number violating processes fall out of thermal equilibrium at the 
corresponding critical temperature. 

A quantitative understanding of the electroweak phase transition is a basic prerequisite 
for the discussion of any model of baryogenesis at the weak scale. This includes reliable 
knowledge of its order and of the strength of the phase transition, if it is of first order. 
Several approaches have been used to investigate the electroweak phase transition. Important 
results have been obtained by use of 3-dimensional effective theory || [7[], e-expansion ||, 
average action || and lattice simulations |T0| , |IT , 12| . In particular, the detailed lattice data 
for physical phase transition parameters from refs. [O, O] permit explicit comparison with 
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results from perturbation theory. 

The present investigation is concerned with the extraction of thermodynamic parame- 
ters of the electroweak phase transition, based on the perturbative calculation of the high 
temperature effective potential, i.e. the free energy of the system. 

Perturbative calculations of the potential to the order g 3 , A 3 / 2 suggest a first order phase 
transition for different Higgs models (g denotes the gauge coupling and A the scalar coupling). 
These calculations, based on the one-loop ring summation, have been carried out in refs. [|1^, 
14"! f° r the Abelian Higgs model and in refs. [|15|, |16j for the standard model. Two-loop 



summation has been done to order g A , A in ref. ||17||, where scalar masses have been neglected 



with respect to gauge boson masses, and by use of another approximation in ref. ITS]. The 
results of ref. 11171 include both the Abelian and the non-Abelian case. 



However, there is a need to extend the work of Arnold and Espinosa to a complete 
g 4 , A 2 -calculation. In the present investigation, assuming formally A ~ g 2 and keeping the 
full dependence on the Higgs field (p, its zero temperature vacuum expectation value v and 
the temperature T, the necessary corrections are calculated |1| p0[ . To obtain information 
about the importance of specific non-Abelian effects and to have a particularly simple model, 
the Abelian case has been considered first [ |L9] . 

Besides the above extension of a conventional approach, a gauge invariant calculation of 
phase transition parameters is presented [21, 331. Using this method a better understanding 



of the physics of the phase transition is obtained. The numerical results of the gauge invari- 
ant approach are similar to conventional Landau gauge calculations, thus supporting their 
reliability. 

Chapter [l] starts with some well known facts about the effective potential, its thermody- 
namic interpretation, its loop expansion and the relevance for the description of first order 
phase transitions. After that the resummation of masses, necessary at high temperature, is 
discussed. Here the emphasis is on a method based on Dyson-Schwinger equations [19|, ^ 
to be used in the sequel. 

In chapter ^ the complete g 4 , A 2 calculation of the finite temperature effective potential 
is performed for the Abelian Higgs model ]nj. This includes a zero temperature renormal- 
ization in the on-shell scheme. The absence of a linear yj-term, explicitly verified at this 
order, is shown to survive to all orders. Latent heat, surface tension and jump of the order 
parameter are calculated as functions of the Higgs mass from different approximations to 
the effective potential. The higher order scalar corrections, added to the previous results, 
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are found to be important at not too small Higgs masses. 

The above analysis is extended to the standard model in chapter || which is based on 



ref. [|20|1 . First, the case of the pure SU(2)-Higgs model, which is much simpler, is discussed 
in detail. The importance of specifically non-Abelian contributions is stressed and the un- 
certainties of perturbation theory are traced back to infrared problems. A qualitatively 
similar situation is found for the complete standard model. The main quantitative change 
is introduced by the large top mass, which reduces the strength of the first order phase 
transition. 

In the last chapter the gauge invariant approach to the phase transition, suggested in 
ref. ||21fl , is described. The phase transition parameters obtained from a one-loop calculation 
are compared with the Landau gauge results, showing good agreement for not too large Higgs 
masses. Several additional questions of the analysis of the electroweak phase transition are 
discussed in the sequel [[22j . They include two-loop resummation problems, the connection of 
gauge invariant and Landau gauge approach, and the use of the Clausius-Clapeyron equation 
in the present context. 

After the discussion of conclusions to be drawn from the above investigation several 
analytic formulae are listed in the appendices. Appendix [A] contains some useful integrals. 
The explicit analytic results for the finite temperature effective potential in the Abelian 
Higgs model, the SU(2)-Higgs model and the standard model are displayed in appendix 
Self energy corrections at zero temperature, necessary for the on-shell renormalization, are 
given in appendix |Cj. 
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Chapter 1 

The effective potential 



1.1 Definition and loop expansion 

The present investigation is concerned with equilibrium thermodynamics, which, as is well 
known, can be completely described as soon as the partition function 



trexp 



-j3(H + [ Jcp) 



(1.1) 



of the system under consideration is given. Here the inverse temperature is denoted by /3 and 
H is the Hamilton operator. The source J is coupled to the field ip, specified later on as the 
Higgs field, which is used as the order parameter for the description of the phase transition. 
In this section it is sufficient to consider the simplest possible case of a (//-model with one 
degree of freedom. The generalization to more complicated field theories is straightforward. 

The fundamental thermodynamic potential per unit volume, W(J), is related to the 
partition function by 

Z = exp(-(3QW) , (1.2) 

where ft denotes the three dimensional volume of the physical system. The temperature 
dependent effective potential is now defined by the transition from the variable J to the 
variable <p, realized by a Legendre transformation: 

V(ip,T) = W(J,T)-Jip , p= dW{ d j T) - (1-3) 



It is straightforward to derive the two identities 

< <p >= if 



dW 1 

<^>= v = -- , -<H > -TS = V(<p,T) (1.4) 



dJ ' ft 

clarifying the physical interpretation of the effective potential as the free energy of the 
system. Here the brackets < ... > symbolize the thermal expectation value of an operator 
and S = — dW I dT is the entropy. 
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The perturbative calculation developed later on is based on the path integral represen- 
tation of the partition function |^3| 



Z = tr exp 



(3{H + J Jtp) = J -D^pcriod. exp J dr J d 3 x(C - Jcp) 



;i.s) 



Here C is the Euclidean Lagrangian and the path integral is taken over all fields periodic 
in time direction. Compactification of Euclidean time results in the replacement of the well 
known loop integral by the integral-sum 

. rp n=+oo „ 

VJ* S (L6) 



n=— oo 



and correspondingly in the introduction of the Euclidean propagator with discrete Matsubara 
frequencies k 

1 ' , k = 2nTn. (1.7) 



k 2 + m 2 k 2 + k 2 + m 2 
This imaginary time formalism, used exclusively in the following, is the simplest formulation 
of finite temperature field theory and perfectly suited for the investigation of the desired 
equilibrium parameters. 

While Z is the sum of all Feynman diagrams, /3QW = — In Z contains the connected 
graphs only. The potential V, obtained from W by means of a Legendre transformation, 
is the momentum independent part of the generating functional of one-particle irreducible 
diagrams. It has to be evaluated at non-vanishing external field, thus requiring the summa- 
tion of infinitely many graphs at each loop order. This problem is solved by the following 
identity 



V(0) = VW*) + i^ft ]»(*> + ml) + {one-particle irred. vacuum diagrams}. (1.8) 

Here V tree denotes the tree-level potential, i.e. the momentum independent part of the free 
Hamiltonian. The curly bracket represents the sum of all one-particle irreducible vacuum 
graphs calculated from the 'shifted' Lagrangian 

C^(ip) = £(0 + ip) — | terms, constant and linear in ip\ . (1.9) 

Correspondingly, m| is the mass term generated after this shift. 



The perturbative methods described in section |L3] are based on the above representation 
of the potential. Note that in the following sections, the shift (p will often be denoted by tp 
for brevity, if no confusion is possible. 
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In the remainder of this section a short derivation of the identity ( |1.8|) shall be given. To 
keep the notation compact, let (3 = ft = 1 during these manipulations. 

Shifting the integration variable in eq. (|1.5|) by some arbitrary but fixed function (p results 



in 



W[J] = -In J Dip exp J C(<p + (p) - {<p + <p) J 
where the Lagrangian after the shift can be written in the form 

C(ip + p) = -V tree {(p) - ipg{(p) + £ v (y?) . 



;i.io) 



n.ii^ 



Note that J is not necessarily constant and W is a functional. The effective action is defined 
by 



T[0] = -W[J] + J JCp 



sw 
17 



;i.i2) 



where arguments have been dropped for brevity in the functional derivative. Introducing 
the notation J = J + g(0) and using the generating functional of the 'shifted' theory 



W*[J] = -\nj Dtpexp J C^(<p) - 



;i.!3) 



the effective action takes the form 



r[$ = -v tree ((p) - w*[J) + (tp — 0)j. 



Observe that 



_ 6W[J] _ 5W[J] _ + 8W*[J] 



;i.!4) 



;i.!5) 



SJ 8J 5J 

Now the shift tp, which has not yet been specified, is set to (p — tp. This gives an explicit 

result for the effective action: 



-V tree (0) - W^[J] at 



5Wv 



o. 



;i.i6) 



Since the effective potential is given by minus the momentum independent part of T, the 
identity ( |1.8| ) follows. 

1.2 Description of a first order phase transition 

The intuitive picture of a first order phase transition, described by the free energy density 
of the system, is based on the double well structure of the free energy as a function of the 
order parameter. In this picture the barrier between the two minima is responsible for the 



8 



Fig. |l.l| Schematic graph of the coarse-grained free 
energy, f(c), and the corresponding true free en- 
ergy, /(c). The dotted sections denote the ana- 
lytic continuation of f(c) into the metastable re- 
gion (from Langer, 



necessity of an activation energy for a transition between the phases, thus rendering the 
phase transition first order. However, a rigorous definition of such a non-convex free energy 
is not trivial. First the phenomenological discussion of Langer shall be briefly described: 
Assume the existence of a coarse-grained free energy, F(c), of the form 

1 



F c) 



2 (dc) 2 + f(c) 



(1.17) 



with some function f(c) characterizing the homogeneous state and an order parameter c of 



the system. The free energy density /, typically of the form given in fig. [O], is well suited 
for the description of the metastable and unstable region. However, even in the region where 
/ has a positive second derivative it is not identical with the true free energy density /. The 
later one is a convex function and connects the two different physical states of the system by a 
straight line. In the metastable, though not in the unstable region, the analytic continuation 
of f(c) does still describe the thermodynamic properties of the one-phase physical system. 

The perturbatively calculated high temperature effective potential is by definition an 
analytic function of the order parameter if. Therefore, assuming its convergence to the true 
free energy in the stable region [T2^, ^7] , it can be naturally interpreted in the above sense as 
the free energy of the metastable states. This, however, does not clarify the interpretation 
in the unstable, non-convex region of the potential. 

Following Langer, the coarse-grained free energy can be calculated by integrating out 
the short wavelength components of the microscopic variable only. This corresponds to 
the introduction of some infrared cutoff characterizing the coarse-graining size. Note that 
the one-loop results for the effective potential are not very sensitive to a small enough 
infrared cutoff. In fact, in the Abelian model no infrared problems are expected at any 
loop-order (see section |2.3|) . Even the non- Abelian two- loop results of section |3.3.1| do not 



change qualitatively when a small infrared cutoff is introduced. This could be taken as a 
justification to interpret the non-convex region of the obtained effective potential in the spirit 
of the coarse-grained free energy. A detailed discussion of the coarse-grained free energy in 
high temperature field theory can be found in ref. . 



At one-loop level such an interpretation is supported by the analysis of ref. |29||. There it 



is shown that in some approximation the effective potential V((p) gives the energy density of 
a homogeneous state with wave functional concentrated on configurations near the classical 
value (p. 

It has to be admitted that the understanding of the perturbatively calculated effective 
potential in the non-convex region does not seem to be satisfactory. This implies some 
doubts about the physical interpretation of quantities calculated from the potential in that 
region. In particular, the surface tension, as it is defined and calculated in chapters and |3] 
is affected by this uncertainty. 

Nevertheless, critical temperature, latent heat and jump of the order parameter are, in 
principle, calculable by standard methods from the effective potential. This is due to the 
fact that they can be obtained from the minima of the free energy, which, following the 
above discussion, are trustworthily described by the potential. 



1.3 Resummation 

The loop expansion of the potential given by eq. (|1.8| ) formally corresponds to an expansion 
in coupling constants. Unfortunately, due to the effectively three dimensional integrals 
arising from contributions with vanishing Matsubara frequencies in the high temperature 
theory terms proportional to (T/m) n do appear. Owing to the Higgs mechanism masses 
are proportional to coupling constants in the relevant theories. Therefore a formal loop 
expansion does not generate the desired expansion in coupling constants. This problem is 
well known and can be solved by resummation. In the following this shall be described in 
detail for some unspecified theory with a generic coupling constant g, to be used as the 
expansion parameter: 

Consider a general Lagrangian with interaction terms generating 3- and 4-vertices pro- 
portional to g 2 and 3-vertices proportional to gk^. Here is a momentum variable, as it 
appears at the vertex in gauge theories. All masses are considered to be of order g. Note 
that this structure is suggested by the standard model Lagrangian, where the square root 
of the scalar coupling a/A, the Yukawa coupling gy, the electroweak gauge couplings (71,(72 
and the strong gauge coupling g s play the role of the generic coupling g. The Abelian Higgs 
model and the SU(2)-Higgs model fit this general structure as well|14], [16|, |30[ . 



Resummation now means that the leading self-energy corrections of order g 2 T 2 have to 
be added to the mass squares, thereby collecting an infinite series of graphs with increasing 
powers of (gT/m). Having realized that in a systematic way, the remaining higher loop cor- 
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rections are always connected with higher orders in the couplings, which ensures a systematic 
expansion of the potential [T(| . 

Working in that spirit, the main task is now to identify all graphs contributing to the 
order g 4 . 

1.3.1 Method based on Dyson-Schwinger equations 

In this subsection a method for the calculation of the effective potential based on Dyson- 
Schwinger equations is presented [[19], 2(J. A similar way of summing the different contribu- 



tions to V for the <£> 4 -theory has been considered in ref . [31] . 

To circumvent the combinatoric problems of resummation it is useful not to calculate 
the potential V(<p) itself but its derivative with respect to the field if, i.e. the sum of all 
one-particle irreducible one-point functions. This 'tadpole' method has been suggested in 
refs. |32|. 



The Dyson-Schwinger equation |33| for the effective potential has the form 

--^(V-V tne ) = A + B = , 



;i.i8) 



where the internal lines represent all particles of the relevant theory and the external lines 
stand for the shifted scalar field. The two different sorts of blobs are full propagator and full 
3-vertex respectively. Therefore the first term has the explicit form 

A = trW(y)rf 2 k , m ^ ■ ( L19 ) 

T k z + mf ree + H(k) 

In general, mass, self-energy and vertex W are matrices and "tr" denotes the sum over the 
suppressed indices. The (^-dependence of the mass m tree , introduced by the shift (|1.9|) , is 
obvious. 



No corrections are needed for the vertex function in the second term of eq. ( 1.18 ). This 
can be verified by the following argument: Any correction to that vertex corresponds to a 
proper three-loop contribution to the potential. Here by proper, the absence of self-energy 
insertions, which would require resummation, is meant. The three-dimensional contribution 
to the three-loop graph, arising when all Matsubara frequencies vanish, results in a g 5 - 
correction. This can be seen most easily by scaling all loop momenta according to k — > kg. 
Similar arguments do also prove that the parts of three-loop diagrams with non-vanishing 
Matsubara frequencies do not contribute 92-dependent terms to the potential up to order g 4 . 
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This way of argumentation is explained in more detail in ref. ||30|| . Appendix A of ref. [|TI 
contains a proof of the sufficiency of self-energy resummation for the order (^-potential, 
which can be generalized to higher orders. 

The self-energy insertions needed in both terms of eq. (|1.18|) can again be obtained from 
a Dyson- Schwinger equation, which, to the required leading order, reads 

-n(A;) = -(II a (A;)+II 6 (fc))= . (1.20) 



It is easily verified that the omitted parts would result in proper three-loop terms in the 
potential which do not contribute to the order g A . In the following the indices 2 and 3 
denote contributions of order g 2 and g 3 respectively. The tadpole part of the self-energy can 
be written as 

n o (A0 = n a2 + n a2 (£0 + n a3 + --- , with n a2 (o) = o. (1.21) 

Here the momentum dependent part n a2 (fc) does only arise in the case of a non-Abelian 
gauge theory. It is introduced by the corresponding projection operator when calculating 
the longitudinal self-energy of the gauge boson (see section |3.1| ). 

The third order term in g does not contribute in the case of nonzero ko- If ko = 0, by 
the above scaling argument k — > kg only its momentum independent part Il a 3 is needed. 

The leading order momentum independent g 2 -part of Tlb{k) will be called n^. All these 
contributions to the self-energy U(k), needed below, are obtained by iterating eq. ( |1.20| ) 
twice. 

Using these definitions and introducing the corrected mass term m 2 , 

m 2 = m 2 tree + II a2 + n fe2 , (1.22) 



equation ( |1.19|) can be written as 

dk 



A = tr 



k 2 + m 2 + n a2 (fc) + n a3 + u b (k) - n 



1/2 



trw(<p)rtdk[ La , * ^ , L , - -J— i n a3 - 



k 2 + m 2 + U a2 (k) k 2 + m 2 k 2 + m? 

Here the second equality is obtained by expanding the integrand in g. Separate treatment 
of the parts with ko = and ko ^ together with the above scaling argument show that all 
contributions of order g A are taken into account. 



12 



Inspection of the last term of eq. (|1.23|) and term B of eq. ( |1.18|) shows that their sum is 



equal to the derivative —dV e /d<p, where — V e represents the sum of all two-loop diagrams 
of the type shown in fig. [L^.a (setting sun diagrams). To see this, notice that 

)A%) = -^m 2 , (1.24) 

a direct consequence of the shift generating both mass terms and 3-vertices from the inter- 
action Lagrangian. This observation simplifies the actual calculation significantly. 
Furthermore, the contribution 

^fVtrWMf^n^ (1.25) 

is treated separately. It is equal to the sum of all terms bilinear in masses coming from 



two-loop diagrams of the type shown in fig. [1.2|.b, i.e. the sum of their three-dimensional 



parts. Note that here and in the setting sun diagrams discussed above, masses resummed to 
leading order have to be used. 



Fig jl.2| Typical two-loop diagrams 



The remaining part 

V R = -f ^trW^fok L + .J +nAk) + ^n S2 ^) (1.26) 



can be easily rewritten as 

I r p cilx, I c elk, 

v R = .j dm H^ WT - ? + _ tr ^ FT - ?( ii - (*) -n B ) . (i.27) 

Therefore the complete potential is given by 



V = V tree +V e + V z + V R . (1.28) 

Denoting by V3 the sum of the tree level potential and the g 3 -order part of Vr and calling 
V4 the fourth order corrections of Vr, 

V 3 + V 4 = V tree + V R , (1.29) 

the following final formula is obtained : 

V = V 3 + V 4 + V e + V z . (1.30) 

This representation, used for the explicit results in the appendix, has the advantage of a 
separation of the third order part from the higher order corrections. 
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1.3.2 Counterterm method 



Another possibility of a systematic resummation at high temperature is based on the in- 



troduction of thermal counterterms |34|. The simplest form of this method is to add and 



subtract a temperature dependent mass term for every field in the Lagrangian. Let 

£ = 4in + £mass + £/ with C mass = - \ m \w\ (1-31) 

i 

be the Euclidean Lagrangian depending on several fields (p^ with kinetic term, mass term 
and interaction term. The leading temperature dependent mass corrections are given by 

-5m t 2 = -IT 2i (0)= , (1.32) 



where only the g 2 T 2 -pieces of the diagrams are considered. The Lagrangian is now rewritten 
according to 

£kin mass J , 



C! mas, = ~ Yl \( m li + <*™i , C'l = £l + j I] 6m2 iV 2 i 

i i 



;i.33) 



The main point is now in the cancellation occurring between the thermal counterterms from 
the new interaction Lagrangian and those contributions from higher loop diagrams which are 
of low order in the couplings and require resummation. For example, the sum of the diagrams 



from fig. [T73| does not contribute to the field dependent part of the effective potential up to 
order g A . Here the dot symbolizes the thermal counterterm. If the thermal counterterm is 
considered as increasing the formal loop order by one, all of these diagrams are of three- 
loop order. Treating the proper three-loop diagrams, i.e. those without self-energy parts, 
separately, it can be shown that the sum of all three-loop diagrams does not contribute to 
the desired g 4 -potential. 



Fig. 1.3 Three-loop diagrams with self-energy parts 



Arguments of that kind can be generalized to higher loop orders. Therefore, to obtain 
the full g 4 - result in this approach, the diagrams shown in fig. |1.4| are sufficient. Here the 
circle symbolizes the one-loop contribution which is proportional to the logarithm of the 
determinant of the propagator. It has to be kept in mind that the notation in this chapter 
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is merely a generic one, so that the lines in the diagrams stand for all the different particles 
of the theory. 

Fig. |1.4| Typical diagrams contributing to the g 4 -potential in the counterterm method. 



In ref. |T7]] a slight modification of this method is used. There, the resummation is applied 
only to the zero Matsubara frequency modes. Thereby many constant terms and terms which 
would cancel each other in the final result are omitted from the beginning. Also, using this 
method, it is not necessary to keep the dependence on the space dimension in the thermal 
counterterms, if dimensional regularization is used. The details of this will not be discussed 
here, because the following calculations are based on the method of section |1.3.1| . However, 
the results of chapters |2| and |] have been checked by use of the counterterm method. 
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Chapter 2 

Abelian Higgs model 



The simplest gauge theory with spontaneous symmetry breaking, the Abelian Higgs model, 
is believed to exhibit the main features of the electroweak phase transition |S5|, . At one- 



loop order the effective potential shows a first order phase transition at high temperature, 
driven by a term cubic in the vector mass in complete analogy to the standard model case 



13| , p~4|| . The study of the Abelian Higgs model might prove useful for the understanding of 
the electroweak phase transition, because this simple model does not suffer from the severe 
infrared problems of the non-Abelian theory. No magnetic mass is expected to arise. This 
does not present a problem due to the absence of the higher-loop non-Abelian contributions, 
which are divergent in the symmetric phase. 

In this chapter a complete calculation up to the order e 4 , A 2 , with gauge coupling e and 



scalar coupling A, is presented, supplying the results of ref. [17] with scalar corrections. An 



analysis of the phase transition parameters shows that these corrections are important if the 
Higgs mass is not too small. 

2.1 Calculation of the potential 

The effective potential is calculated by expanding it in the coupling constants, as described 
in section |1.3.1| . Consider the Euclidean Lagrangian 

C = --F^F^ - \D^\ 2 + z/|$| 2 - A|$| 4 , (2.1) 

with 

D^d^ + ieA^ and $ = -j=(<p + tp + i X ). (2.2) 
Using the tree level vacuum expectation value v the Higgs mass term is rewritten as v = Xv 2 



and counted as order A. Then the identification with the the generic coupling g of section O 
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reads g ~ e ~ vA. At this point it is necessary to describe in more detail the resummation 
procedure for the vector particle Hl4 , ^3] : 



The choice of the Landau gauge, which will be used throughout this investigation, is 
justified by the absence of vector - Goldstone boson mixing together with the ip- independence 
of the gauge. Two-loop calculations for the three-dimensional model do also show that the 



convergence of the perturbation series is best in Landau gauge [[37]]. The bare propagator 
can be given in the form 

D {k) = ¥ ^(P L + P r ), (2.3) 

where 

3 / kk\ k k 

Pt^u = J2 S »i y S ij ~ -j^J S jt> alld P L^ = <W ~ P T»v (2.4) 

Here the full covariance has been broken down to SO (3), which is the symmetry left after 



the specification of a rest frame intrinsic in the formalism of thermal field theory. In ref. [[14 
the full propagator is shown to have the form 

D(k) = ±D 0{ k)m)D ( k )r = j?T - ? L__ ft + _J__P T , (2 .5) 

thus permitting a simple derivation of the mass corrections from the self-energy U(k). 
Note that Pl and Pt are orthogonal projection operators. 

Now the calculation of the different contributions to the potential, listed in eq. ( |1 . 28| ) , 
will be described in some detail. This is possible due to the extreme simplicity of the model. 
The explicit results are found in appendix [BTT[ Throughout this investigation dimensional 
regularization is used. This means that the naively three dimensional parts of the integral- 
sums are evaluated in n — 1 = 3 — 2e dimensions. This section together with the relevant 
parts of the appendix describe the calculation of the MS-potential which will be improved 



by the zero-temperature renormalization of section ^2 



The first step is the calculation of the leading order temperature corrections to the masses. 



Here basically the results of ref. |L4J have to be supplemented with the e-dependent parts 
which produce finite contributions to the potential due to one-loop divergences. The results, 
which can be obtained using the integrals of appendix [A.l| , are given in appendix [B.l . 



Consider the contribution Vr first. Due to the particle content of the theory, given by 
Higgs particle, Goldstone boson, transverse and longitudinal vector boson, the first term of 
eq. ( |1.27|) reads 

Vbx = g {/ dm 2 v I(m v ) + J dm\l{m x ) + (2 - 2e) J dm 2 T I{m T ) + J dm 2 L I{m L )} , (2.6) 
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with the standard temperature integral 



4 dk wh*- (2 - 7) 



Km, 

k z + mr 

Similarly, the second term of eq. ( |1.27| ) is given by 

V R2 = -^(2 - 2e)IL j2 , T J(m T ) - ^n fc2 , L I(m L ) . (2.8) 

Note, that this type of corrections does only exist for the vector particles and that U a2 (k) = 
in the case of the Abelian model. The sum of both terms Vr = Vr\ + Vr 2 enters the 
contribution V4 which can be found in appendix jB.l , 



The two-loop contributions of V z are of vector-scalar and of scalar-scalar type. Figure 



271] shows the setting sun diagrams of the Abelian Higgs model, contributing to V e . The 
labelling follows the standard model case. As usual, dashed and wavy lines represent scalar 
and vector propagators respectively. When solving the integrals, it is useful to write the 
propagator as a sum of a covariant part and a longitudinal correction 

DW = B » w + (FT<-FT^) Pi ' (2 ' 9) 

This is possible, because the transverse mass receives no leading order temperature correc- 
tion, due to gauge invariance. 

Fig. |2.1| Setting sun diagrams in the Abelian Higgs model 

A straightforward calculation of all the terms described above results in the explicit formulae 
given in the appendix. 

The necessary MS-counterterms are generated by a multiplicative renormalization of 
couplings, mass term and scalar field 

X b = Z x \ , e 2 b = Z e 2e 2 , v h = Z v v , §\ = Z^ 2 , (2.10) 

where the index h denotes the bare parameters. No renormalization of the vector field is 
required because it does not enter the final formula for the potential. The counterterms for 
the potential are now easily obtained by expressing the bare parameters, which enter the 
order g 2 , A part of the potential, through the MS-quantities, defined by eq. ( |2.10| ). Note that 



the corrections 5Z = Z — 1 are defined to have no finite part. They are given explicitly in 
appendix |B.1| together with the correction to the potential which they generate. 
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2.2 Renormalization at T = 



To get rid of the arbitrary scale j2 the potential is rewritten in terms of physical parameters 
defined at zero temperature. Such parameters are the Higgs and vector masses and the 
vacuum expectation value of the physical Higgs field (Pp hya = Z~ 2 (fl- Note that here, in 
contrast to the MS-definitions of the previous section, the counterterms 5Z do have finite 
parts, to be specified below. 

Returning for the moment to Minkowski space, the usual on-shell definitions of field 
renormalization factor, Higgs mass and vector mass read 



dq 



ml + Re n v (mj) = m 2 ^ phys , m 2 + Re n(m 2 ) = m 2 phys . 
They are completed by the specification of the vacuum expectation value v 

dV 



= 0, (2.12) 

tPphys=V 



d(Pphy 

where the Landau gauge zero temperature effective potential to the order e 4 , is given by 



u b 2 A b 4 m \ 1 1 3 /i 2 
2 4 b4ii z \ e 2 m 



V = V tree + V 1 = -^i + ^t-^ + +l n . 



(2.13) 



G I tt- \ r 2 m 2 j 647T 2 [ e 6 m- 



The one-loop self-energies of Higgs particle and vector field in Landau gauge are given in 
appendix |C1| . Note that in this approach, no one-particle reducible tadpole contributions 
need to be considered because they vanish due to eq. (|2.12|) . 
Requiring the validity of the tree-level relations 

m l, P hys = 2 ^ 2 , m \h V s = eV , v = Xv 2 (2.14) 

at one-loop, the above equations permit a straightforward derivation of the necessary coun- 
terterms: 
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5Z^=Re%{ml) , 5Z X = -Re%(mD + ±- (-V{(v) - Ren^m* 

(2.15) 

5Z V = — (-V{(v) - Ren^(mJ)) , 5Z e2 = -Ren;(mJ) - ^Ren(m 2 ). 

The correction to the potential can now be obtained by renormalization of its leading order 
part, i.e. by inserting these counterterms into eq. ( |B.3| ). This results in a potential, explicitly 
independent of the renormalization scale p. However, in view of the principal features of the 
potential considered here, the numerical effect of the performed finite renormalization is not 
very important (see fig. |2.4Q. 



2.3 Absence of a linear term 

In the early days of the perturbative treatment of the phase transition in the U(l)- and 
SU(2)-Higgs models, the possibility of a linear yj-term in the order-g 3 -potential has been 
discussed. It has however been realized that no such term is present, if the resummation is 



performed correctly |14|, [38|, |39], [40|| . Of course, in the present calculation the cancellation 
of linear mass terms to order g 3 is reproduced. These terms are not displayed in the final 
formulae in the appendix to make them more compact. 

However, the presence of linear terms in higher orders has not been completely clarified 
before. Their cancellation is claimed in ref. [[40] in a variational approach and in ref. by 
some gauge invariance argument which is not further specified. Therefore it is interesting to 
observe that in the present result the linear <£>-terms cancel to order g 4 : Consider the explicit 
formulae of appendix [B.l| at temperatures above the barrier temperature Tj, defined by the 
vanishing of the (^-independent scalar mass term. Expanding these expressions in f at the 
point if = 0, linear terms are found in V a and V z , which cancel each other exactly. This 
feature, formulated more precisely as 

dV 

lim — — = to all orders in e and A , 

<p^o dip 

can be shown to survive to all orders of small couplings resummed perturbation theory. The 
proof is based on an identity following from global U(l)-symmetry: 

1 dV / \ 
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Obviously it suffices to demonstrate the finiteness of the self energy n x (g 2 = 0) in the limit 
ip — > 0. Due to the positive temperature masses of x an d <p singularities can only arise from 
the transverse gauge boson propagator above the barrier temperature. Therefore diagrams 
of the kind shown in fig. |2.2j have to be investigated. Here the wavy lines symbolize leading 
order resummed vector propagators and the blobs are full vertices without internal vector 
lines, meaning the sum of all diagrams built from scalar propagators with the correct number 
of external vector lines and possibly one or two x-lines. Notice that the vector resummation 
affects only the longitudinal modes, and is therefore irrelevant for the discussion of small-y? 
singularities. 

Fig. |2.2| Higher loop self-energy corrections to m 2 

Consider first the full vertices with external vector lines only denoted by 

rS^C*!,..., (2.17) 

below. Since in the contributing diagrams all propagating particles are massive scalars, the 
vertices F are analytic in fcj. 

Above the barrier temperature scalar masses are analytic in ip 2 . In addition, explicit 
y?-factors appearing at the vertices of some diagrams are always paired, due to the structure 
of the unbroken theory, which has no vertices with an odd number of scalar lines. This shows 
that r is also analytic in p 2 . 

Furthermore the three-dimensional part of T with vanishing external Matsubara frequen- 
cies satisfies 

rS...^(fci, • . • , k 2n , p> = 0) ~ \h\ . . . \k 2n \ (2.18) 

for small \ki\ and k® = , a/3 . . . \w E {1, 2, 3} . 

This follows from a gauge covariance argument, completely analogous to the zero temperature 
case. Having established these properties of T, the small-y? behaviour of the x-seh-energy 
can be derived as follows: 

Consider the most dangerous lowest power of p stemming from the maximal infrared 
divergence, which is obtained by setting k° = for all transverse vector propagators. It can 
be calculated by scaling the loop momenta according to k — > kip. The above discussion of 
the pure vector vertices V shows that after this scaling they can be counted as order ip 2 at 
least. In the case where no scalar line connecting two external x _ hnes exists (see the first 
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diagram of fig. |2.2j) , the power counting in (p proceeds as follows: a factor p 3 for each of 
the L loops, p~ 2 for each of the I internal vector lines, p 2 for each of the V — 2 full vector 
vertices and an explicit factor p for each of the two vertices with a x-line. Together this 
gives the minimal over all power of p 

n v = 3L-2I + 2(V -2) + 2 (2.19) 

for a diagram with V vertices (compare the argumentation in appendix A of ref. fl6fl). 

If there is at least one scalar line connecting the two external x-lines (compare the second 
diagram of fig. [2.2|) , V — 1 full vector vertices contribute factors p 2 . In this case however the 
last term +2 does not exist because a vertex with two x-lines need not have an explicit 



(^-factor. Therefore eq. ( |2.19| ) is valid in the second case as well and consequently in general. 



Now the well-known formula V + L — I — 1 immediately gives 

n v = L>0, (2.20) 

or equivalently: There is no divergence for (p — > 0. 

If some of the vector propagators have non-zero Matsubara frequencies, the vertices 
connected by those "heavy" lines may be formally fused. Now repetition of the above 
argument leads again to the desired result, thus completing the proof. 

This nice feature of the Abelian model strongly supports the hope for a reliable pertur- 
bation series in the symmetric phase. Unfortunately, due to the 3- and 4-vector vertices of 
the non- Abelian theory, the above argument does not apply there. 

2.4 Numerical results and discussion 

Fig j2.3| Different approximations of the effective potential plotted at their respective critical 
temperatures at mHi ggs = 38 GeV (the e 4 , A-potential is a result of ref. [ |I~7|j ) 



Fig. 2.4| Dependence of the surface tension calculated from the different potentials on the 



zero temperature Higgs mass 

Fig j2.5| Influence of the e 4 <^ 4 -correction on the surface tension at small Higgs masses 

For the numerical investigation of the phase transition standard model values are chosen 
for vacuum expectation value v and vector boson mass mw at zero temperature (see sec- 
tion |3.3.1p. Of course the analogy has to be used with caution, because even at one-loop 
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level, the three vector degrees of freedom of the SU(2) increase the strength of the phase 
transition considerably if compared to the U(l) case. Nevertheless, assuming the standard 
model situation, in the following analysis the Higgs mass will be considered as an unknown 
parameter. 

As is well known, the potential suggesting a first order phase transition is generated 
at order g 3 , A 3 / 2 by the combination of quadratic, cubic and quartic yj-terms, where the 
coefficient of the quadratic term is of the form T 2 — T 2 and therefore strongly depends on 
the temperature (see V3 in appendix |B.1|) . Figure |2.3| shows the different approximations 



to the potential at their respective critical temperatures and mniggs = 38 GeV, all of them 
suggesting a first order phase transition, but of quite different strength. Calculations of 
order e 3 , A 3 / 2 |TB|, [CJJ], of order e 4 , A JH] and the present e 4 , A 2 -calculation are compared, 



showing a dramatic decrease of the barrier height in both higher order results. 

Obviously, reliability of perturbation theory has to be questioned already at this small 



Higgs mass. A more detailed picture can be obtained by considering the surface tension [|41 

s = £ + d<p^2V(^T c ), (2.21) 

which may be seen as a measure of the strength of the phase transition. The reliability of 
the present calculation of this quantity is of course doubtful, because it relies essentially 
on the non-convex region of the potential, where the physical interpretation is still unclear. 
Nevertheless s can be used conveniently to discuss the properties of the potential as a function 
of the Higgs mass. The results are shown in fig. |2.4j . 

Fig. |2.6| Higgs mass dependence of the latent heat AQ of the phase transition 



Fig. |2.7| Position of the second minimum ip + at the phase transition in units of the critical 



temperature T c 

The reasons for the difference between the e 3 , A 3//2 and the e 4 , A 2 -results are twofold. Con- 
sider the region of small Higgs masses first. Here the e 4 </9 4 -term, being a large correction to 
the tree-level term A<^ 4 /4, is mainly responsible for the decrease of the barrier height. This 
is illustrated in fig. |2.5| , where in addition to the consistent third and fourth order results a 
potential containing the third order terms together with the e 4 (/? 4 -correction (see V4 in ap- 
pendix |B.1| ) is investigated. The effect is not removed by zero temperature renormalization, 
although the temperature independent constant c\ forms the largest part of the coefficient of 
e 4 ip 4 . Notice that this large constant arises from the expansion of the temperature dependent 
part of the one-loop integral I(m) fi2"|| , and is therefore absent at T = 0. 
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As it can be seen from fig. |2Jj|, two-loop contributions are not too important at small 
Higgs mass and perturbation theory appears to be in a relatively good shape. 

For Higgs masses above ~ 30 GeV another higher order effect becomes more important: 
The three-dimensional pieces of two-loop temperature integrals generate terms of the form 
(p 2 \n(m + m^^) and the like (see V a and Vj, in appendix |B.1|) , which, in spite of their 
numerical smallness, influence the potential significantly. This can be understood by recalling 
that at the critical temperature the leading order <£> 2 -terms essentially cancel and that a in- 
dependence in a coefficient of ip 2 can not be absorbed in a correction of T c . These logarithmic 
terms with a positive sign decrease the barrier height, which is clear from the shape of 
the function x 2 ln(x + const.). In summary, the most infrared sensitive contributions of 
the high temperature field theory introduce large corrections and prevent the reliability of 
perturbation theory. 

The effect of the above logarithmic terms is overestimated by the approximation used in 
ref. \ V7\ . There, scalar masses are counted as order A 1//2 and neglected systematically against 
vector masses, resulting in contributions of the type ip 2 In m. These terms, lacking the scalar 
mass cutoff in the logarithm, destroy the first order phase transition for Higgs masses above 



~ 40 GeV. The plot at 38 GeV of fig. |0| does already show the arising pathology of the 
e 4 , A-potential, which can be somewhat eased but not cured by addition of the A 3 / 2 terms. 

The latent heat of the phase transition is another interesting quantity to be calculated 
from the effective potential : 

&Q = Tjj j y{<p + ,T)\ Ta , (2.22) 

where ip + is the position of the asymmetric minimum of the potential V, normalized to zero 
at the origin. This relation follows easily from the definition of AQ together with the formula 
relating entropy and free energy: 

AQ = TAS = T (5 symmetric Shroken) i S = 7^; . (2.23) 

Figure [2.6| shows that the change of the latent heat, introduced by the g 4 , A 2 -corrections, is 
very large, but not as dramatic as for the surface tension. The phase transition appears to 
be much weaker first order at two-loop. 

The vacuum expectation value in the broken phase at T c , shown in fig. |2.7| , does not 
reflect the dramatic change of the surface tension, introduced by higher order corrections. 
This quantity, the jump of the order parameter, seems to be the most reliable characteristic 
of the phase transition, accessible in perturbation theory. 
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Note that the critical temperature is, here as well as in the standard model case, very 
close to the uncorrected barrier temperature, defined by the vanishing of the quadratic term 
in eq. ( |B.7|) . Therefore a graphic representation does not seem advisable. 

Depending on the considered quantity and the standards to be chosen the phase transition 
can be regarded as understood in principle for Higgs masses below 30... 50 GeV. For larger 
masses it is likely to be much weaker first order, although even such a qualitative description 
is not really well founded due to the infrared problems. The next sections will show that 
reliability at physically relevant large Higgs masses remains a problem in the non-Abelian 
theory as well. 
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Chapter 3 
Standard model 



In the following standard model discussion the technical parts will be confined to the new 
features not present in the Abelian model. While the non-Abelian character of the theory 
changes the numerical effect of two-loop corrections qualitatively, fermions and the additional 
U(l)-symmetry are less important. Therefore it suggests itself to concentrate on the much 
simpler SU(2)-Higgs model for a detailed analysis of the two-loop effective potential. 



This chapter is based on ref. [20 



3.1 Calculation of the potential 

To fix the notation the essential parts of the standard model Lagrangian are given: 

t-'Higgs t~-gauge £ fermion ^Yukawa- (^•^) 

Defining the covariant derivative as 

Dp = d, + igJ^Bp + %yPK; (3.2) 

the gauge part and fermionic part for nj fermion families are unambiguous. The Higgs 
contribution reads 

W = -P,»l" + ^I»-W • where « = -L ( _ £ + ^ j (g 3) 

denotes the Higgs doublet and the shift has been applied according to tp\ — > (p + tp\. 
All fermions except the top quark are considered to be massless. The resulting Yukawa 
Lagrangian is 

Cvukawa = -9YqL§t R , q L = ( *^ , $ = ir 2 $*. (3.4) 
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Following section |1.3j the formal power counting rule 

9i ~ 92 ~ 9y ~ A 1/2 (3.5) 

is used. 

The resummation of the scalar and vector degrees of freedom proceeds in analogy to 
sections |1.3.1| and [2.1| . Explicit formulae are found in appendix |B.2| . The vector resummation 
is complicated by the mixing of the W3- and 5-fields, characterized by the well known 
Weinberg angle 9 in the transverse part and by a different, temperature dependent angle 
9 in the longitudinal part. To implement this in the evaluation of Feynman diagrams the 
notation of a transverse and longitudinal propagator is introduced: 

T k 2 + m 2 T ' L k 2 + m 2 L L ^ 
Now the WVWVpropagator is conveniently rewritten as 

D W 3 W 3 = D Z cog 2 q + D J gin 2 q + ^ cog 2 q + ^7 gin 2 £ (37) 



= cos 2 + Dl sin 2 + (cos 2 ^ - cos 2 9) [(D T + Dfj - (D T + ,D2) 

Here the masses are specified by the indices Z and 7. From this the 5-Z?-propagator is 
obtained by exchanging sine and cosine. Similarly, the -B-H^-propagator reads 

Several comments are in order concerning the calculation of Vr. Calculating the first term 
of eq. (|1.27|) the e-dependent part of ~rn^ L J rm\ L is needed (compare eq. (|2.6|)). It is obtained 
most easily from the last line of eq. (|B.18j ). For the second part of Vr, corresponding to the 
second term of eq. (|1.27|) , the self energy contribution U a2 (k), introduced in eq. (|1 . 21 



is 



required. It is nonzero for the iy-field only and can be found in appendix |B.2| together with 
the contributions of type Hb 2 . These terms are displayed most conveniently in the original 
.B-WVbasis. 

V z contains two-loop contributions of vector- vector, vector-scalar and scalar-scalar type. 
The setting sun diagrams are shown in fig. |3.1|, where the labelling follows ref. |T7j|, to make 



comparison more easy. 

V e = V a + V b + Vi + Vj +V m + V p . (3.9) 

Note that ghost contributions have been included in V m , and that the purely scalar diagram 
V p has not been considered in ref. [171 . The calculation of all these terms is long but 
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straightforward. All the integrals with complicated covariant structure can be reduced to 
the basic types given in appendix |A.1| , as described in ref. |L7 . 



Dropping the appropriate terms of V the lower order g A , A-result, as it is given by Arnold 



and Espinosa in ref. |17| , can be derived. This is also valid for the Abelian model discussed 
in the previous chapter. Another calculation for the SU(2)-Higgs model appeared in ref. 0, 
where the complete g 4 , A 2 result is derived from the effective three-dimensional theory to- 
gether with some corrections from ref. jTT). The analytic results of the actual loop-calculation 
are found to be in agreement with the present analysis. However, due to a specific way to 
apply the renormalization group method the numerical outcome of ref. [M] differs from the 
results to be presented here. 



Fig. 3^ Setting sun diagrams for the standard model 



Note that there are linear <£>-terms of fourth order in the couplings present in V a , V m and 
V z . As in the Abelian model, these terms cancel each other, thus ensuring the relation 
lim^o dV/ dip = for all allowed temperatures. This cancellation is essentially the same 
effect which leads to a vanishing third order transverse gauge boson mass term in the sym- 
metric phase ]16j, as can be seen in the contributions of diagrams fig. 6.o, 6.q, 6.t and 6.u 
of ref. ||16| . 



3.2 Renormalization at T = 

In the Abelian Higgs model the performed zero temperature renormalization has not affected 
the calculated phase transition parameters significantly. Nevertheless, it is not possible to 
profit from that experience by just setting /2 = 1/(3 in the standard model case. The reason 
for that is the large negative gy</9 4 -term which dominates over the tree level quartic term. 
This leads to an MS-potential unbounded from below for small Higgs mass. 

The zero temperature renormalization is performed in the on-shell scheme, as described 



in ref. [ |3] . Slightly modifying the procedure of section |2.2| , the physical parameters chosen 



are Higgs mass, top quark mass, W- and Z-boson masses and the fine structure coupling a, 



defined in the Thompson limit P4| . The physical masses and the wave function renormaliza- 
tion of the Higgs field are defined in analogy to eq. ( [2.11| ). A multiplicative renormalization of 
the coupling constants, the tree level Higgs mass square —v and the Higgs field is performed. 



The required one-loop corrections are of course well known in Feynman-'tHooft gauge [T3| . 
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However, here they are needed in Landau gauge. There is no problem with the correction to 
the electric charge 5e = e& — e, which is gauge independent [fj5| . This can be easily checked 



explicitly using the results of ref. ||46|| , where the gauge dependence of several self-energy and 
vertex corrections has been calculated. Therefore in the present calculation the formula for 
Se from [43| is used. The logarithmic terms with the five light quark masses are treated in 



the way described in ref. |47|], with data from ref. [[4§], resulting in the vacuum polarization 
contribution 



Ren^(Ml) = Ren£(M|) - < a >) = -0.0282 ± 0.0009 



(3.10) 



The dependence of the one-loop self energy corrections on the gauge parameters has been 
calculated in ref. [pi for gauge bosons. Therefore the corrections in Landau gauge, needed 
here, can be taken from [|6[ ff9|| . The self energy corrections for the physical Higgs boson 
and the top quark can be easily calculated in Landau gauge. The results are displayed in 
appendix 



Defining v as in the Abelian model by eq. (|2.12|) , no one-particle reducible tadpole dia- 
grams need to be added to the self energies. The zero temperature one-loop effective potential 
Vtree + V\ , required for that, is easily obtained from eq. (|2.13| ) by including the additional 
degrees of freedom. The renormalized couplings and the renormalized Higgs mass square are 
defined by 

c=^- , C g l = sg 2 = e , A = — , M H = 2v , (3.11) 



Mz oivi w 

where s = sin 8\y and c = cos6w Now the counterterms follow easily from the definitions 
of the physical parameters: 



A7 J e n w 

0Z„2 = Z FT H 7T- i 

31 e Ml Ml ' 



' w 



5Z„2 



e s 2 M z 



c 2 n 



w 



s 2 M 2 z 



M 2 W 



M 2 H 



+ 



92V( 



c 2 U w 

s 2 M 2 v ' 



2 ' 



2M W M- H 



5Z V 



SZ. 



Suit 1 



UY 



m t 



-- 2 K 



(3.12) 



M 2 H 



+ 



2M W M 2 H 



Here n and n' stand for the real parts of the self-energies and their derivatives at the on-shell 
point. 

The renormalized, /2-independent potential can now be obtained by applying eqs. (|3.12| ) 
and the field renormalization dZ^ = H' to the leading order contribution, given by the first 
line of eq. (|Fmp . 
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Clearly, the resulting formula for the potential is too long to be given explicitly. However, 
it seems worthwhile to give the numerically most important parts of the corrections to enable 
a simplified usage of the analytic result in the appendix. As it has already been mentioned, 
the main contributions come from the g^-corrections to parameters of order A (see also [l7|): 

3gy , mt . 3g Y v 2 
— !r In -^7- , ov — 

?7T /i 16n A 



5A = ^ ln — > 5Z/= 17-^- ( 3 - 13 ) 



Introducing these corrections in all terms in the potential contributing to order A and using 
standard model tree level relations to calculate the couplings one obtains a result which 
is 'partially renormalized at zero temperature'. The corresponding correction to the MS- 
potential reads 

sv = t(- 6u+ w 6X ) + 5 ^- (3 - 14) 



As it will be seen later (section |3.3.2| ), the numerical effect of this drastic simplification is 
not too severe. 

3.3 Numerical results and discussion 

Fig. |3.2| The surface tensions calculated from the different potentials as functions of the zero 
temperature Higgs mass 



Fig.3.3 Influence of the most infrared sensitive contributions on the surface tension as a 



function of the Higgs mass 

Fig. |3.4| Higgs mass dependence of the latent heat AQ of the phase transition 



3.3.1 SU(2)-Higgs model 

To obtain an understanding of the qualitative effects of higher order corrections the pure 
SU(2)-Higgs model is studied first. In this subsection the additional U(l)-symmetry and the 
effect of fermions are neglected. A discussion of this simplified version is also useful in view 
of lattice investigations, which deal with the pure SU(2)-Higgs model presently and probably 
also in the near future. 

The relevant potential can be easily derived from the formulae given in the appendix 
|B.2| by performing the limit gi,gy - and setting the number of families rif to zero. This 
results in the potential of appendix B.3 , which is the basis for the numerical investigation of 
the present subsection. 
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Standard model values for W-mass and vacuum expectation value v are used, unless 
stated otherwise : My/ = 80.22 GeV and v = 251.78 GeV. The parameter p, of dimensional 
regularization is set to T = 1/(3. This can be justified by the small dependence on the 
renormalization procedure. The differences between the results obtained in this scheme and 
in a scheme with on-shell T = renormalization are relatively small. This phenomenon is 
observed in the Abelian Higgs model as well. 



In analogy to the investigation performed in section 2.4 the potentials from calculations 



to order g 3 , A 3 / 2 |L5], [L6[], to order g 4 , A JT7| and to order g 4 , A 2 p(J are compared. All suggest 



a first order phase transition in a wide Higgs mass range. The form of the potentials at 
the critical temperature is the standard one (see fig. |2.3| ) and will not be shown here again. 
However, comparing the barrier height in the different approximations the picture differs 
significantly from the Abelian case in the region of large Higgs masses: Both g 4 , A- and 
g 4 , A 2 -potential suggest a much stronger first order phase transition than the lowest order 
result. The pathological behaviour of the e 4 , A-potential of the Abelian model does not arise. 

In the following, the surface tension (see eq. (|2.21|) ) will be used to illustrate the features 
of the potentials at different Higgs masses. Figure |3.2| shows the similarity of the situations 
in the Abelian and the non-Abelian case in the small Higgs mass region. The g 4 (p 4 -term 
from the one-loop vector contribution is responsible for the strong decrease of the barrier 
height in the higher order results. 

Fig.|3.5| Position of the second minimum (p + in units of the critical Temperature T c 



Fig j3.6| Surface tensions from the third and fourth order potentials for a model with 
M\y = 20 GeV in the Higgs mass range where the values differ by a factor of 2 at most 

At larger Higgs masses, approximately above 40 GeV, the infrared two-loop contributions 
become more important. Their effect is however quite different from the Abelian case. 
Compare the results of order g 3 , A 3 / 2 and g 4 ,\ 2 first. The increase in the strength of the 
phase transition, studied already in ref. [|18|], can be traced back to the infrared features 



of a non-Abelian gauge theory. The crucial contribution is the one coming from the non- 
Abelian setting sun diagram (fig. |3.1| .m). It produces contributions to the potential of type 
ip 2 ln(f3mw) with a negative sign. Recall, that the logarithmic terms from diagrams |3.1| .a and 
[3~T|.b (or ^I|.a, |2j].b), discussed in the Abelian case, have a positive sign. Notice also, that 
the new terms are non-analytic at </? = 0. The reason why these kind of corrections affect 



the form of the potential strongly has already been discussed in section 2.4. To demonstrate 
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Fig.|3.7| The surface tensions obtained from the standard model effective potentials as func- 



tions of the zero temperature Higgs mass with mt op = 170 GeV 

Fig. |3.8| Position of the second minimum ip> + in units of the critical temperature T c in the 
case of the standard model with m top = 170 GeV 

its importance the ip 2 ln(/?m H /)-term of V m has been deleted by hand. The corresponding 
surface tension is shown in fig. |3.3| (long-dashed line). 

The comparison of the potentials to order g 4 , A jl7| and g 4 , A 2 shows a picture very similar 
to the Abelian Higgs model. Neglecting scalar masses in the logarithms results in spurious 
(p 2 ln(/3mw)-terms, which reduce the surface tension. Another important contribution, less 
significant in the Abelian case, is the one proportional to g 2 (mi + 3rri2)mwL from V z . This 



term comes from scalar- vector diagrams of type of fig. |1.2| .b and it was neglected in ref . [17] . 
On the relevant scale (<p < T) it produces a very steep behaviour of the potential, again 
increasing the surface tension. The observed difference between the result of ref. |I7j and 



the complete g 4 , A 2 calculation presented here is mostly due to these two effects, together 
with the well known influence of the cubic scalar mass contributions from V3. However, in 
sharp contrast to the Abelian case, both the g 4 , A- and the g 4 , A 2 -calculation suggest a much 
stronger first order phase transition than the lowest order result. 

Another interesting effect of higher order A-corrections is the complete breakdown of 
the phase transition at a Higgs mass of about 100 GeV, where the surface tension is very 
large. In this region the above mentioned term, proportional to g 2 {m\ + 3m2)mwL, becomes 
important. For a temperature close to the uncorrected barrier temperature TJ,, at which the 
scalar masses vanish for ip = 0, it produces an almost linear behaviour in the small (p region. 
This results in a potential for which at T = Tb the asymmetric minimum is not a global 
minimum but only a local one. Note that Tf, is the lowest temperature accessible in this 
calculation. In other words, the temperature region in which the phase transition occurs can 
not be described by the given method, due to infrared problems. 

In order to illustrate the possible effects of the unknown infrared behaviour of the trans- 
verse vector propagator, the dependence of the surface tension on the magnetic mass can 
be studied. A magnetic mass of order g 2 T/37r is motivated by the solution of gap equations 
TBJ and supported by the numerical investigation of gauge invariant gap equations with 



resummed vertices ||50|| . Following the approach of ref. |T6| the transverse vector mass takes 
the form 



"4=(^v+(?4Y. (3.i5) 



2 J V 3vr/3 
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where 7 is some unknown parameter. One can introduce this redefined transverse mass 
in the most influential infrared contributions, i.e. in the mfp- and in the tp 2 In /Jm^-terms. 



Figure |373| shows the results obtained for 7 = 0, 2 and 4, supporting the qualitative behaviour 



found in ref. ||16|| . The main difference is due to the fact that the higher order result suggests 
a stronger first order phase transition, thus for a given mHi ggs a larger magnetic mass is 
necessary to change the phase transition to second order. 

A complete fourth order calculation of the surface tension has to include the wave function 
correction term Z v (ip 2 ,T) calculated in ref. |p8|| . Using the results of ref. |28j the surface 
tension 

a = r + d<pyj2(l + Z v (<p*,T))V(<p,T c ) , (3.16) 

J 

has been determined for Higgs masses between 25-95 GeV. The numerical effect of this Z- 
factor is very small, only 1% — 4%. This is due to the smallness of the potential in the only 
region where Z^ is significant, i.e. at small ip. 

The latent heat AQ (see eq. (|2.22 ))is plotted in fig. |3]4] as a function of the Higgs mass. In 
contrast to the Abelian model, here the higher order results show an almost linear increase of 
the latent heat with the Higgs mass. This somewhat surprising behaviour can be understood 
by observing that for those potentials neither the position of the degenerate minimum nor the 



height of the barrier change significantly with increasing Higgs mass (see fig. |3.2| ). Therefore 
dVj dT does not change dramatically over a wide Higgs mass range. On the other hand the 
critical temperature increases with growing m H i ggs - 



As in the Abelian case the quantity cp + /T c , shown in fig. |3.5| as a function of the Higgs 
mass, is least affected by higher order corrections. 

Now the question arises whether a good convergence of the perturbation series, which 
can not be claimed in the whole range of A for a realistic gauge coupling g = 0.64, could 
be present in the region of much smaller gauge coupling constants. This seems indeed to be 
the case, as can be seen in fig. |3.6| , where the surface tensions of order g 3 , A 3//2 and g 4 , A 2 are 
plotted for a model with a vector mass of 20 GeV, i.e. g = 0.16. In the used Higgs mass 
range the two results for s differ by a factor of two at most. The relative size of this range, 
i.e. the ratio of the minimal and maximal values of the Higgs mass, is 4, which is twice as 
large as the corresponding range for the model with My/ = 80 GeV. 

The phase transition parameters calculated in this section can be compared to new lattice 
data available at a low Higgs mass point (mni ggs ~ 18 GeV) and a high point (mni ggs ~ 
49 GeV) 0, |i~2| . These data include critical temperature, jump of the order parameter 
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latent heat and surface tension at both Higgs mass values. The g 4 , A 2 -results are in 
good quantitative agreement with lattice data (explicit numbers will appear soon fl22f). This 



is highly non-trivial, since quantities change by large factors between low and high point. 
An exception is formed by the surface tension at the high point, which is larger by a factor 
of ~ 2.5 in perturbation theory. All other parameters agree between the perturbative and 
the lattice calculation with deviations compatible with the observed scaling violation on the 
lattice and the uncertainties of perturbation theory. Note, that the simulations have been 
performed with a maximum number of three lattice points in time-like direction, thus the 
results may change somewhat on larger lattices. To obtain the good agreement quoted above 
the renormalization effect of the vector boson on the very light Higgs particle has to be taken 
into account at the low point 



3.3.2 Complete standard model 

In the case of the full standard model the qualitative behaviour of the potential is essentially 
the same as for the SU(2)-Higgs model. The main difference is a decrease of the surface 
tension by a factor ~ 4. This can be traced back to the large top quark mass which 
influences the potential by lowering the barrier temperature (see |B.2ip . The additional 



U(l)-symmetry and the light fermions are less important. Also the characteristic points of 
the surface tension plot of fig. are shifted to higher values of the Higgs mass. Figure 
37r| shows the surface tension as a function of mHi ggs - The complete breakdown of the g 4 , A 2 



calculation, observed at m H i ggs ~ 100 GeV for the pure SU(2) case, occurs at mni ggs ~ 
200 GeV in the full model. These quantitative differences do not change the qualitative 
features of the potential, thus the discussion given in the previous section does also apply 
to the standard model. The difference between the fully renormalized potential and the 
partially renormalized potential (see eqs. ( |3.13| ),( |3.14|) ) is not too severe in view of the huge 
uncertainties still present in the perturbative approach. Again, the position of the second 
minimum at the critical temperature given in fig. |3.8| , does not depend as strongly on the 
order of the calculation as the height of the barrier. Unfortunately, the region mHi ggs ~ 
40 GeV, in which the reliability of the perturbative approach is the best and <£>+/T c w 1, is 
well below the experimental Higgs mass bound. 
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Chapter 4 

Gauge invariant approach 



In this chapter the gauge invariant treatment of the electroweak phase transition is presented 
following ref. [^TJ. The main new point is the introduction of a gauge invariant source term 

PI 



Z = exp(-pttW) = trexp 



-P(H+ ( J&$>) 

Jv 



(4-1) 



and the subsequent definition of an effective potential by the Legendre transformation 

V(v,T) = W(J,T)-J°- , | = ^™, (4.2) 

which is performed perturbatively. The new variable a/2 =< > describes the thermal 
expectation value of a gauge invariant quantity. 

This approach is perfectly well suited for comparison with lattice investigations, which 
usually proceed without gauge fixing and consider the expectation value of the operator 
Note also, that effective actions for composite operators have been defined previously, 



for example in ref. [p2|1 , and that a formulation based on a field linearly related to (ft has 



been given in ref. [53fl. 

4.1 One- loop calculation of W(J) for the 
SU(2)-Higgs model 

The analysis is restricted to the pure SU(2)-Higgs model, which is sufficient to illustrate the 
main features of this method at the one-loop level. Note that a comparison with the Abelian 
model makes no sense at one-loop level since qualitative differences start to appear only at 
the next loop order. Using the Lagrangian £ from section |3.1| , where the appropriate limits 
have been taken (see section p.3.1|) , the thermodynamic potential W can be calculated from 
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the path integral 

exp(-/ftW) = J D$D&DW^exp ^ dr J d 3 x (£ - J$ f $) 



(4.3) 



To evaluate this integral the extremum of the static, ^-dependent part of the action with 
source term, characterized by V tree , has to be determined: 



C J = C - J$t$ = -V tree + 



V, 



tree 



+ J + a T 2 )&$ + \(&<f>Y 



(4.4) 
(4.5) 



Here, following the philosophy described in section |1.3.2| , a thermal counterterm for the 
scalar field has been introduced. Its leading and next to leading order contributions read 



«o = «oi + "02 



1, 3 2 

-A H q 2 

2 16 y 



3g 3 



(4.6) 



167T V 6 

The g 3 -term, generated by the one-loop self-energy contribution of the longitudinal vector 



boson (see ref. |16 |, eqs. (41), (42)), is taken into account for reasons to be explained below. 

Two different regimes have to be distinguished for the calculation of W(J). Consider 
the case with J < v — aoT 2 first. In this situation Vtree develops a non-trivial minimum at 
ip i = 0(J) (see eq. ( |3.3|) for the conventions), defined by 



\ip 2 = M 2 = v - J - a T 2 . 



Performing the shift tpi — > <p + (pi a tree-level contribution 

W tree (J) = ±(-v + J)£ 2 + ^ 4 



(4.7) 



which is independent of the resummation, is generated. The one-loop contribution in Lorentz 
gauge (cf. ref. [Q) with gauge fixing parameter a reads 



W X {J) = ^dk 



9 \n(k 2 + m 2 ) + \n(k 2 + mj) + 3 ln(P + k 2 m 2 2 + am 2 m 2 2 ) - 6 In k 



+ 



d?k 



3 ln(A; 2 + m\) - 3 ln(fc 2 + m 2 



(4.9) 



2 J (2tt) 3 /5 

Here Higgs mass, Goldstone mass, vector mass and longitudinal vector mass are given by 



m \ = 2M 2 , mi = , m 2 



- 2 -2 

:9 <P , 



1 f> 

2 2^2 i 2t-i2 

™ L = -g v +-gT . 



(4.10) 
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Note that only the zero Matsubara frequency mode of the longitudinal vector degree of 
freedom has been resummed. This simplifies the calculation in Lorentz gauge and does not 
change the result up to order g 3 extracted from the one-loop formula. 

If, on the other hand, J > v — aoT 2 , no shift of the field is necessary and W(J) has no 
tree-level contribution. The one-loop term can be easily obtained from eq. ( f4.9|) by setting 

m x =m 2 = M with M 2 = -v + J + a T 2 (4.11) 



and replacing transverse and longitudinal vector masses by zero and m^o = \j5/QgT respec- 
tively. 

In both the symmetric and the broken phase eq. fl4.9|) gives explicitly gauge independent 
results since the product is always zero. Working in i?^-gauge the same answer is 

obtained. The integrals are easily performed using appendix |A.1| and terms of fourth and 
higher order in the masses are neglected together with constant terms common to both 
phases. 

The one-loop thermodynamic potential W( J) in broken and symmetric phase is given by 
W(J) = W b {J)Q{u - J - a T 2 ) + W s (J)e(-u + J + a T 2 ) , 



where 



W b (J) 



T 2 
—M 2 
6 



12vr 



—M 

AX 

r2\3/2 



(4.12) 



(4.13) 



(2M 2 Y /2 + Qm 6 + 3ml 



W S (J) = ^-M 2 
o 



T 

127 



4M 3 + 3m| , 



(4.14) 



and M 2 is given by eqs. ( ^4.7| ) and ( |4.11| ) respectively. 



4.2 Gauge invariant effective potential 

The one-loop gauge invariant effective potential can now be obtained by a perturbative 
Legendre transformation according to eq. ( |4.2| ) . In the usual approach, based on the order 
parameter if, the perturbatively defined effective potential is the sum of the one-particle 
irreducible vacuum graphs of the shifted theory. However, no such interpretation is known 
for the gauge invariant potential V(a). Therefore, after writing W(J) and J as a sum of 
contributions of increasing order 



W{J) = W (J) + WtiJ) + W 2 {J) + 



J = J + J l + J 2 + 



(4.15) 
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the definition (|4.2| ) has to be implemented order by order in perturbation theory |54 
The first terms of the resulting potential are given by 

2 



V(a) 



W Q (J 



5* 



Wi(J„ 



+ 



win 1 r^( J ) 



•A) , 



<9J 2 



Jo 



+ 



where Jq is a function of a defined by the leading order relation 



dW (J) 



d.J 



a 

Jo 2 



(4.16) 



(4.17) 



Consider the phase with broken symmetry first. Here W(J) has to be decomposed into 
a leading term Wo (order g 2 , A) and a next to leading term W\ (order g 3 , A 3//2 ), defined by 
the first and second line of eq. (|4.13|) respectively. 

In the symmetric phase, where no tree-level term exists, the Legendre transformation can 
be performed exactly at one-loop order. Neglecting constant terms common to both phases 
the resulting potential is given by 

T 2 



V(a) = V s (a')e(a') + V b {a')Q{-a') 



a = a 



3 ' 



(4.18) 



where 



Via') 



\{a 0l T 2 -u)a' + \a' 2 



(4.19) 



T 



12tt 



{2Xa'f 12 + 6 Q$V) 3/2 + 3 (\<?o> + \g 2 T 2 



3/2 



Vs {a') = \a T 2 -u)a'-- < ^- — (-g 2 T 2X " 1 
y ' 2 V ; 6 T 2 Arr W 



(4.20) 



Here the new variable a' has been introduced to separate the shift of the field variable 
generated by the interaction from the basic, model independent thermal expectation value 
T 2 /3 p5] . Note the similarity of V& to the order-g 3 , A 3 / 2 Landau gauge result, eq. ( |B.31| ). 
The new values of the scalar masses, namely zero for the Goldstone boson mass and 2Acr' 
for the Higgs mass square, form the main difference of and V3. 

Fig j4.1| Gauge invariant effective potential at critical temperature, mHi ggs =70GeV 



Fig. |4.2| Effective potential in Landau gauge at critical temperature, mHi ggs =70 GeV 



The gauge invariant potential defines a critical temperature at which its two minima, 
one in the symmetric and one in the broken phase, are degenerate. This is illustrated in fig. 
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4~l| . For comparison, the Landau gauge potential to order g 3 , A 3 / 2 is plotted at its critical 



temperature in fig. |4.2| . The interaction induced shift a' seems to play a role similar to the 
squared field expectation value (p 2 of the usual approach (compare the discussion in section 



The main difference of the gauge invariant and the Landau gauge potential lies in the 
description of the symmetric phase. Choosing the expectation value of the basic field oper- 
ator as order parameter the symmetric phase corresponds to exactly one point, where this 
parameter must vanish. In contrast to this, the gauge invariant coupling of the source per- 
mits a description of the symmetric phase by a non-trivial minimum of the free energy as a 
function of the field square expectation value. This expectation value can be made smaller 
than its value at the minimum by turning the source, i.e. the mass term, on. Therefore the 
gauge invariant potential exhibits a very steep behaviour left from the symmetric minimum, 
but no definite smallest value of a'. 

Fig. |4.3| Latent heat as a function of the zero temperature Higgs mass, calculated from the 
potential in Landau gauge and from the gauge invariant potential 



Fig. |4.4| Shift of the Higgs field as a function of the Higgs mass, calculated from the potential 



in Landau gauge and from the gauge invariant potential 

Note, that the gauge invariant potential is continuous together with its first derivative at 
the matching point a' = 0. This is due to the difference of the first terms of eqs. Q4.19Q and 
(|4.20| ) which compensates the contribution of the last term of eq. (|4.19| ) to the derivative at 



a' = 0. The additional resummation up to order g performed in section [4. 1| is responsible 
for that matching. However, the same result for V(a) can also be obtained with leading 
order resummation only. In that case the contribution 

3o 2 T 2 1 

-j^Mm^ = -a 02 T 2 a' (4.21) 



from the scalar-vector diagram of the type of fig. O.b appears in the symmetric phase, thus 
rescuing the above matching of the first derivatives. 

At small Higgs masses the symmetric minimum of the gauge invariant potential is not 
important numerically. Therefore in this region the form of the potential at the critical 
temperature and the derived phase transition parameters are very close to the Landau gauge 
results. However, already at mHi ggs = 70 GeV the symmetric minimum leads to an increase 



of the barrier height by a factor of ~ 2 (compare figs. and |4.2|) . This effect becomes 



39 



even more important for larger Higgs masses. The latent heat plotted in fig. [O] does also 
show a stronger first order phase transition at large Higgs mass values suggested by the new 
approach. Note that a similar behaviour of the latent heat as a function of the Higgs mass 
has been obtained from the two-loop results in Landau gauge (see fig. |3.4|) . The jump of 
the order parameter shown in fig. |4.4| is not seriously affected by the new approach. This 
justifies the conjecture of section |3.3.1| that this parameter is reliably determined by the 



Landau gauge calculation. 

Due to the separate treatment of the symmetric and the broken phase in the gauge 
invariant approach, it is not clear at present how the surface tension can be calculated. 
However, due to the strong increase of the barrier height at large Higgs masses in the new 
approach significant changes with respect to the Landau gauge results are expected. 

4.3 Problems at higher orders 

It is, in principle, straightforward to extend the calculation of the previous sections to the 
order g 4 ,X 2 , as it has been done in the conventional Landau gauge calculation of chapters ^| 
and [| Nevertheless, the concrete realization of this project is hampered by some problems 
to be discussed in the sequel. 

The first step is to supply the leading order results for W(J), given by eqs. (|4.13|) 



and ( 4. 14]) , with higher order corrections. It is advantageous to calculate in Landau gauge, 



although of course the final result should be gauge independent by definition. The necessary 
corrections include the next term of the high temperature expansion of the one-loop integrals 
and the leading contributions of the two-loop graphs. The later ones, which form the main 
part of the calculation, can be obtained using the formulae of appendix |B.3| and changing 
the masses appropriately. Note that in the broken phase, one-particle reducible two-loop 
graphs have to be considered. 

To make W(J) finite, the parameters of the Lagrangian are renormalized multiplicatively 
in the conventional way, using e.g. the SU(2)-limit of eqs. flB.20| ). Note that the multiplica- 
tive mass renormalization has to be applied to both — v and J. However, a finite result is 
obtained only after the additional counterterm 

W c . t . = -L-(J - vf (4.22) 

has been added to W(J) in both the symmetric and the broken phase calculation (compare 
ref. [0). It is claimed that a possible finite part of this counterterm, common to both phases, 
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does not affect the physical parameters derived from the calculation. Firstly, constant and 
linear term in J do only shift the Legendre transformed function along vertical and horizontal 
axis respectively. Applying such a shift to both phases does not change the description of 
the phase transition. Secondly, the quadratic term in J does not affect the minimum of the 
Legendre transformed function, since at the minimum J = and the quadratic correction 



vanishes together with its first derivative [53]. To see this more explicitly introduce the 
corrected function W(J), defined by 

W(J) = W(J) + CJ 2 . (4.23) 

The new position of the minimum of the potential is given by 



o min dW 



2 dJ 

and analogously 



dW 



j=o dJ 



= , (4.24) 
j=o 2 



V mm = W(0) = W(0) = V mm . (4.25) 

Performing the Legendre transformation perturbatively a higher order dependence of the 
result on the finite part of W c ,t. is nevertheless present. 

Now the obvious way to proceed is to calculate V(a) using eq. (|4.16 ). In the broken 



phase W2 is given by the order g 4 , A 2 corrections and the first three terms of the perturbative 
expansion of V can be calculated. In the symmetric phase, where no tree-level term exists, 
these two-loop corrections define W\ and only the first two terms of eq. (4.16|) are to be used. 



It has been checked that the complete potential is continuous at the matching point 
a' — 0. However, its first derivative is logarithmically divergent at this point. Since all 
physical information is extracted from the minima of the potential this pathology should, in 
principle, have no importance. 

The convincing form of the potential illustrated in fig. |4J] has been obtained due to 
the additional resummation of the scalar mass. More precisely, the scalar mass has been 
resummed to the same order to which the whole calculation was performed. Therefore it ap- 
pears necessary to extend the resummation to the order g 4 , A 2 in this section. Unfortunately, 
this can not be done in a straightforward manner, because the two-loop scalar self-energy is 
divergent at zero momentum. Another possible way to organize the higher order resumma- 
tion would be to demand the continuity of the first derivative of the potential. This fails due 
to the same infrared divergences. However, it has been seen numerically that such a higher 
order resummation can change the physical picture of the phase transition significantly. 
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It is not known at present, how the above problems, which are essentially connected with 
the infrared divergences of a massless three-dimensional field theory, can be resolved. 



4.4 On the relation to the conventional effective 
potential 

The gauge invariant thermodynamic potential W(J), introduced in the beginning of this 
chapter, has a simple relation to the effective potential, as it can be calculated in any fixed 
gauge |56j : 



exp(-[3nW) = J DQD&DW^exp J dx{C - J$ f $) 

= 2tt 2 J <p 3 dtp J D$D& DW„ 5 (ip-^J^dx <p 

/ dx(C - J&$) 



(4.26) 



i=2 
„3 



2tt 2 J </dy exp [-pnV((p, J)} 



Here cpi are the real components of <3> (see eq. (|3.3|) ) and V(ip, J) is the Landau gauge effective 
potential of a theory with mass square —v + J. In the infinite volume limit only the absolute 
minimum ip m in of V contributes to the (^-integral in eq. ( |4.26| ), resulting in the relation 



W(J) = V(<p mtn , J) . 



(4.27) 



The thermodynamic potential W(J), calculated by this method, is continuous, but its first 
derivative has a discontinuity at the critical temperature and J = 0, generated by the jump 
of the absolute minimum ip m i n from zero to the non-trivial, symmetry-breaking minimum 
<p + . This discontinuity corresponds to the shift of the field square expectation value a 



2 



dW 



dJ 



dW 



j=o. 



dJ 



j=o 4 



d_ 

dJ 



V(ip+,J)-V(0,J) 



(4.28) 



Considering only the leading terms of the Landau gauge potential V, the approximate result 
is 

" 2 Aa = dJ I 2 { ~ V + J + amT2) ^+ + ' ' " 1 w 2^+ ■ (4 ' 29) 
This relation justifies the way in which the Landau gauge and the gauge invariant calculations 
have been compared in section [4.2| , in particular in figures |4.1| , [4.2| and [4.4| . 

The conventional definition of the latent heat from W(J) gives immediately, via eq. (|4.27|) , 
the formula used in the Landau gauge investigation (eq. |2.22| )). 
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Note however, that the gauge invariant W(J) has been calculated in a completely different 
manner in section [4.1| . There, it has only played an intermediate role in the calculation 
of V(a). The approach was based on the identification of the minimum of the tree-level 
potential and a subsequent loop expansion around this minimum. 

In contrast to the paragraph above, a different approach has been taken in this section. 
After shifting the scalar field, the effective potential V((p, J) which includes loop corrections, 
is determined. Its minimum then gives the thermodynamic potential W(J). 

The way in which the scalar masses have been resummed in the Landau gauge calculations 
of chapters [| and ^ is not unique. If only the minima of the potential are to be considered, 
as suggested by eq. ( [4.27D , it may be advantageous to resum differently in the symmetric 



and in the broken phase. In particular such a method can take into account the fact that 
the Goldstone boson mass has to be zero at the broken minimum. One-loop calculations, 



based on this idea, reproduce the numerical results of section |4.2| exactly. A corresponding 



two-loop investigation is under way f22 | 



The above method, based on the derivatives of W(J,T) near its non-analytic point, has 
its disadvantage. No information can be extracted concerning the metastable and unstable 
states, thus disabling the usual derivative expansion approach to the surface tension. 

4.5 Clausius-Clapeyron equation 

In classical thermodynamics there is a well known relation between the latent heat of a phase 
transition and the change of the molar volume, the Clausius-Clapeyron equation [|57|| . In 
complete analogy a relation between latent heat and jump of the order parameter a can be 
written down. 

Describing the state of the system by the gauge invariantly coupled source J and the 
temperature the phase transition curve is given by a function J cr it.(T). Since W(J,T) is 
continuous in the J-T-plane, the total derivative of W along the phase curve 

dW dW dJ crit dW , Ann . 
= _| (4.30) 

dT 8J dT or y 1 

is equal in both phases. The partial derivatives of W can be easily related to the field square 
expectation value a and the energy density E: 

dW 1 dW 1 

= 5 , , W -jiW-E). (4.31) 
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Now, evaluating eq. (|4.30Q in both phases and equating the results gives 

1 dJ^^A^ 

2 dT T y ' 

where AQ = E symmetric - E broken is the latent heat. 

Neglecting effects of higher loops und renormalization dJ cr i t ./dT can be easily evaluated 
by dimensional arguments. Since the Higgs mass term — v + J is the only dimensionful 
parameter of the SU(2)-Higgs model the dependence of the critical temperature on J is 
given by 



T = const. yV- J cr it. ■ (4.33) 

From this 

T dJcrU. = ^ (4 34) 



dT J=0 

and therefore the jump of the order parameter and the latent heat are related by 

AQ = -im Higgs Aa , (4-35) 

where mHi ggs is the zero temperature Higgs mass. 

This relation has been checked against the available data from one- and two-loop cal- 
culations of the effective potential. The Landau gauge results are in very good agreement 
with eq. ( |4.35| ) at small Higgs masses, with deviations increasing up to ~ 25% at m H igg S =80 



GeV. It is interesting to observe that the gauge invariant one-loop calculation of section H72 



satisfies the above relation exactly. This can be proven analytically, by writing the gauge 
invariant potential in the form 

V{a) = --a + T 4 V(a/T 2 ) (4.36) 

with some dimensionless function V, which has no implicit temperature dependence. From 
this general form of V and the definitions of Act and AQ the relation fl4.35|) can be easily 
extracted. 

Note also, that the lattice results of ref. verify eq. (|4.35| ) quite well. For L t =3 
lattices the deviation amounts to no more than 1...2 standard deviations, defined by the 
statistical error of the simulation. Here L t denotes the lattice size in time like direction. 
An improvement of the validity of eq. ( 4.35|) is observed when going from L t = 2 to L t = 3 
lattices. 
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Conclusions 



The thermodynamic parameters of the electroweak phase transition have been analysed in 
a perturbative approach based on the high temperature effective potential. A complete 
calculation of the g A , A 2 -potential has been performed for the Abelian Higgs model, the 
SU(2)-Higgs model and the standard model. The Abelian calculation has been shown to 
have no infrared problems in the systematic coupling constant expansion in the symmetric 
phase. This does not hold in the general case, where the typical non-Abelian diagrams 
become important in higher orders. This is seen explicitly at the order g 4 , where, opposite 
to the Abelian case, logarithmic mass terms increase the strength of the first order phase 
transition dramatically. 

Critical temperature, latent heat, surface tension and jump of the order parameter 
have been calculated for different Higgs masses and different approximations to the poten- 
tial. The reliability of the perturbative expansion is clearly worsening with increasing Higgs 
mass. In particular, in the non-Abelian case, any information on an infrared cutoff, e.g. a 
magnetic mass value, would increase the accuracy of the calculation drastically. However, 
infrared problems are also connected with scalar masses, which become small at the critical 
temperature. 

Comparing the complete standard model with the pure SU(2)-Higgs model no qualitative 
change is found. The main difference is a decrease of the strength of the phase transition 
due to the large top quark mass. 

Newly available lattice data at mmggs ~ 18 GeV and at mmggs ~ 49 GeV [[□], [12| are 
in good quantitative agreement with the perturbative two-loop results. This may be seen 
as a justification to take the calculated parameters more seriously at larger Higgs mass 
values, where convergence of the perturbation series is bad and no lattice data is available. 
Perturbative results suggest a weak first order phase transition at realistic Higgs mass values. 



However, a change to an analytic crossover is also a possibility [p0| . Unfortunately, the 
region where reliable predictions are available and where <p ~ T, an important condition for 
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baryogenesis, is at mmggs ~ 40 GeV and therefore well below the experimental bound. The 
present results seem to disfavour scenarios of electroweak baryogenesis relying on a strong 
first order phase transition within the minimal standard model. 

The gauge invariant description, elaborated in the last chapter, allows a better physical 
understanding of the thermodynamics of the phase transition. In this approach, coupling 
the source term in a gauge invariant manner, a more direct access to physical quantities 
is possible. In particular, the symmetric phase is described by a non-trivial minimum of 
the potential as well. Since the numerical results obtained in the gauge invariant approach 
at one-loop are similar to the Landau gauge results, the latter ones are strongly supported 
by the new, conceptually more satisfactory treatment. Applying the Clausius-Clapeyron 
equation to the electroweak phase transition a simple relation between latent heat and jump 
of the order parameter has been derived in the above context. Being in good agreement 
with perturbative as well as with lattice data, it improves confidence in the correctness of 
the treatment of the phase transition. 

In agreement with recent results of other methods the performed investigation predicts 
a first order electroweak phase transition of decreasing strength when the Higgs mass is 
increasing up to mmggs ~ 70 GeV. At larger Higgs mass values the calculation is strongly 
affected by infrared problems. If, in spite of this difficulty, perturbation theory without 
explicit infrared cutoff is taken seriously, the phase transition is found to remain weakly first 
order. 



46 



Acknowledgements 



I am most grateful to my advisor, Prof. W. Buchmiiller, who suggested this investigation, 
for continuous support and encouragement. During the intensive and fruitful collaboration, 
I had many invaluable opportunities to learn from him. I would like to thank Dr. Z. Fodor, 
with whom I collaborated during large parts of this investigation, for sharing his ideas with 
me, for many highly instructive discussions and, in particular, for patiently explaining the 
lattice results to me. At different stages of the work presented in this thesis, I profited from 
discussions with Dr. D. Bodeker, J. Hein, Dr. T. Helbig, Dr. B. Kniehl, Dr. H.-G. Kohrs and 
Dr. O. Philipsen. I thank them for their interest in my work. Finally, I would like to thank 
all my colleagues and friends at DESY for the enjoyable working atmosphere. 



47 



Appendix A 
Some integrals 



A.l Integrals of thermal field theory 



In this section the basic temperature integrals |17|, |34], [42], |58(], needed throughout the calcu- 



lation, are listed. The notation follows the general formalism introduced in section 1.1 . At 



one- loop level, the fundamental integral is [17, 42] 



I(m) 



i 



dk 



2c 



k 2 + m? (2vr) n - 1 /5 



+ rnr 



(A.1) 



1 



12P 2 



l + a e ) 



m in" fl , ._, ._, 3 



where 

d = § + 2hi47r-27£^ 5.4076. (A.2) 

Here terms of higher order in m/3 have been neglected because in the present calculation 
they do not contribute to the g 4 -potential. The renormalization scale of the MS-formalism 
is defined by 

ln/i 2 = ln/i 2 + ln4vr - 7 B , (A.3) 
where 7^ is Euler's constant. The coefficient 

fi 2 P 2 12 A/ , 



te = 2 TB + ln^ -C'(2) 

4 TT Z 

of e from the leading contribution 



[1 + eO 



(A.4) 



(A.5) 



" up 2 

will cancel in the final formula for the potential, similarly to the temperature dependent 



divergences W7 . 
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The fermionic analogue of this integral is 
J ' (m) = Tf 



(A.6) 

1 m 2 fl , „ 3 



1 + e(t e - In 4)) - — — - + In jl'p~ A + - - Cl + 2 In 4 



24/3 2 v v ' " 16tt 2 Ve ^ ^ 2 
Its leading, mass independent, part will be denoted by Ij^. 

Another bosonic one-loop integral, appearing due to the non-covariant structure of the 
vector propagator, is given by 

ff w -^ffb) - (4 + £ ) n - £ G + ln ^ 2 + 1 - c • (A - 7) 

At two-loop level the basic problem is the calculation of the scalar setting sun diagram 
(see fig. |1.2| .a). This has first been done in ref. [34] and extended to the case of three different 



masses in ref. [T7|. Neglecting terms of order m {3 the integral reads 



H(m 1 ,m 2 ,m 3 ) 



M dk dg 



J J {k 2 + m 2 )(q 2 + rn|)((A; + q) 2 + 77t|) 



(A. 



3 re . 1 (l -2n2 A i P( m l + m 2 + ™s) 

1% + r . 9/39 m/i p - 4 In c 2 



167r 2 e " 64tt 2 /3 2 
where the constant 

c 2 = 2 - ^ + In ~ - 1^ « 3.3025 (A.9) 

has first been obtained analytically in ref. |58| . 

In the standard model calculation the fermionic analogue of this two-loop integral, i.e. the 
corresponding integral where two of the propagators have fermionic Matsubara frequencies 



k = (2n + l)vrT, is needed. However, as has been shown in ref. [IV] , this integral vanishes 
in leading order in dimensional regularization when e — > 0: 

xfxf dkdq 1 2 2 

f/f(k 2 + m 2 ){q 2 + m 2 )((k + q) 2 + m 2 ) = J 2 ^ P >' ( } 

A. 2 One- and two-point functions 

For the convenience of the reader the usual one- and two-point functions, necessary for the 
renormalization at T = 0, are displayed. Working in n = 4 — 2e dimensions only the pole 
and the finite parts are shown. 

a i 2n ^ 2e f d n k im 2 fl 1 , P 2 \ / \ -i -i \ 

Aim 2 ) = / — = - + 1 + In At , A. 11) 

v ; {2n) n J k 2 -m 2 + ie 16n 2 \e m 2 J ' v ; 
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d n k 



B(q\ ml ml) ^ J ^ _ ^ + + ^ _ ^ + . g) 



16tt 2 



1 n i A* 
- + 2 + In 



where 



a; = 



777.Q + m 2 — q 2 
2m mi 



and the function / is defined by 

V^ 2 — 1 (— arccosh(— x) + 27r) : 



/(*) = 



— vl — x 2 arccos x 
\J x 2 — 1 arccoshrr 



x < 
-1 < z < 
1 < x 
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Appendix B 

Explicit formulae for the effective 
potential 



B.l Abelian Higgs model 



The leading order resummed masses of Goldstone boson, Higgs particle and longitudinal 
vector boson read 

e 2 

Y10 l 



m 2 = Ay? 2 - v + (4A + (3 - 2e)e 2 )/| = \ip 2 - v + A \ +J e + 0( 



m% = m 2 x + 2\ip 2 , (B.l) 

ml = eV + 4e 2 (l-e)/| = eV + ^ + 0(e) . 

The transverse vector mass m? = m = e<p remains uncorrected. For the calculation of Vr 
the self-energy parts of the type n 62 (see eq. ( |1.20|) ) are needed. They appear only for the 
longitudinal and transverse vector boson: 

IW = (2 - Ae)e 2 T p , n b2jT = -2e 2 /^ . (B.2) 

The counterterms rendering the potential finite are generated from the leading order contri- 
bution (see the first two terms of eq. (p.7| ) below): 
,2 r 

>{8Z V + 5Z^) + r p {4A(5Z A + 5Z V 2) + (3 - 2e)e 2 {5Z e 2 + 5Z^)} 



5V 



2 



+^\5Z X + 25Z V 2), (B.3) 



where 



,'3e 4 \ e 2 

5Z X = ( — - 6e 2 + 10A j p , 5Z e 2 = —p , 5Z U = (4A - 3e 2 )p , bZ^ = 3e 2 p , (B.4) 
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with the pole in e denoted by 



1 



lb7T 2 e 

The different contributions to the potential, listed below, have to be summed according to 
eq. ( p. .30|) , where 

V e = V a + V b + V p (B.6) 



is given by the diagrams in fig. |2.1| . Linear mass terms, poles in e and terms proportional 
to t e (see eq. (|A.4Q ), which cancel systematically in the final result, are not shown and the 
limit e —>■ has already been performed. A finite contribution from 5V has been added to 
V4. The following terms form the complete analytic MS-potential: 



1 f 1 , 1 



P 2 V3 4 



A 4 



1 



127T/3 



m 3 , + mr, + 2m 3 + m| 



(B.7) 



32tT 2 /? 2 



Jj 1. . Jfl2 1 1 1\ , , 
m V 3 _ 12 ^ _ 6 Cl + 4 C2 + 4/ -m ("V + m x) 



(B.8) 



+m ip m x + ^(w 2 + m x)(~ 41n ^ + hi/z 2 /3 2 - c 2 ) - ^( m y _ m 2 £ ) 2 ln(w ¥3 + m y 



— {m ip -m x im v 
m 



+ m x ) + — jm 4 -2(m 2 + m 2 )m 2 + (m 2 -m 2 ) 2 } ln(m v + m x + m) 



Un 2 f3 2 



5m 2 (~c 2 - In j + f8m 2 - 4m 2 + — | j ln(2m + m^) 



(B.9) 



2 771 

(m 2 — m 2 ) 2 ln(m + m p ) + 4m 2 ln(2m£ + m v ) H j lnm^ + 2mm^ + m 2 



16vr 2 /3 2 



In — - c 2 - ln{mj(m v + 2m x )} 



(B.10) 



eV 
64n 2 (3 2 



l A -e 2 (-+ln// 2 /3 2 - Cl 
3 U8 PP 



(B.ll) 
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64tt 2 



m 4 , + + 3m 4 



ln/Z 2 /? 2 + ~- Cl ) -2/;r 



32vr 2 /3 2 



-A(m 2 + m 2 ) + Xm ip m x + e 2 (m L + 2m)(m v + 



(B.12) 



B.2 Standard model 

The leading order resummed masses of Goldstone boson and Higgs particle read 

: \(f 2 - U + 



V _ v + 12^ + ^ + *^ + + ° (e) 



(B.13) 



m 2 = m 2 . + 2\{p 2 . 

While the transverse vector boson masses and the fermion mass remain uncorrected to 
leading order 



1 la 1 

m w = -92V , rn z = m w / cos9 w , m f = ~^9y^ 



(B.14) 



the longitudinal SU(2)xU(l) mass matrix receives temperature corrections in the diagonal 



elements 15 



5m% L 



9l 



91 



(10- 186)71-871/(1 -e)I 



92 (h _l 1 



| + j n/ j+0(e) 



(B.15) 



(2-2e)II-fn f (l-e)P fp 
They result in longitudinal masses defined by 



91 a , 5, 



/ m 2 w + 5m^ t -\9\gw 2 
\ -\9i92V 2 m% + <Jm| L 



/ cos 9 sin 6* \ / m 2 ZL 



\ - sin 9 cos / V m 



7 L 



cos 9 — sin 



sin 9 cos 



. (B.16) 



In the following the short hand notations 



s = sin 9w , c = cos 9w , s = sin 9 , c = cos 9 



(B.17) 
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are used. From eq. ( |B.16|) several helpful identities can be derived: 



m ZL m ^L 
2 i 2 



m 



c 2 5m 2 BL + s 2 5m 



z 

m 2 7 + bvnLrj + 5m 2 



2 

WL 



+ 5m Rr 5m 



BL UI "'WL 



WL 



U BL 



4 i 4 

m ZL + m. 



m 2 7 + Sm 2 VL + 5m 



BL 



2m| 



c 2 8m 2 BL + s 2 5m 2 VL 



(B.18) 



For the calculation of Vr self-energy parts of the type II a2 (/c) and 11^ are needed: 

kl 



yrscalars 
U b2 



tjWW 
Ll a2,L 



(k) = 2g 2 F p (l-e)$, 



-nWW 
ii 62,T 



tjBB 
1L b2,T 



-$2(6-46)1, 
-9\l% 



-nWW 
f3 i LL b2,L 



ttBB 
LL b2,L 



9l 



91 



2(3 - 8e)l£ - 8n/(l - e)!^ 
(1 - 2e)/f - f n/(l - 6)1 



fP 



(B.19) 



In complete analogy with the Abelian case (see eq. ( B.3f )) the counterterms are generated 



from the leading order potential given by the first line of eq. (|B.21|) . The multiplicative 
renormalization realizing the MS-scheme is described by 



<5 'Z<fP = {^gl + ^9i - 3#y) p, SZ g 2 



9! 



43 4 



SZ X = \ (-30* + ^ + |^ + ^g* + 6\g 2 Y - \\gl - \\g 2 2 + 12A 2 ) p , (B.20) 



5Z, 



9Y 



17 2 9 2 



<? 2 - 4^ P 



5Z„ 



3c/y - - -^g\ + 6A ) /; . 
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As in the appendix |B.1| the final result is presented without linear mass terms, poles and 
t e -terms, which would cancel in the sum. The following contributions form the complete 
standard model effective potential to the order g 4 , A 2 : 



1 f 1 , 3 2 1 2 1 2 

-» + ^{- 2 X + T6 92 + -^ + 4 9y 



A 



-<P 



(B.21) 



127T/3 



mi 3 + 3 m 2 3 + 4 m w 3 + 2 m W L + 2 m z 3 + m ZL 3 + m lL 3 



V a 



g\ 



32n 2 (3 2 



711 



W 



1 1 A L P 1 , 2„2 1 1 1 

2 - — + — In In u 2 (3 2 - -c x + -c 2 + - 

c 2 2c 4 / 1 3 12 PM 6 4 4 



(B.22) 
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+ ^) | ( m i + 3m l) In ^ + In /i 2 /3 2 - c^j + 2m 2 (mi + m 2 ) 

-4s 2 m 2 In (2m 2 ) - — ^ — f 1 + (mi - m 2 ) 2 {mi + m 2 ) + 2m 2 m z s 2 
2i7iw V 2c/ 



— 2 (m 2 + m 2 ) mly + (m 2 — m 2 ) j In (mi + m 2 + mjy) 



2m 



2 1 m W 



|m| — 2 (m 2 + m 2 ) m| + (m 2 — m 2 ,) | In (m x + m 2 + m z ) 



+ 



(ic 2 ~~ ~~ ln ^ 2m2 + + \ ( m w ~~ 4m2 ) ln(2m 2 + m w ) 



V h 



9l 



64vr 2 /3 2 



c 4 + 4 + 4c 2 + 4 - 10 ) m 2 w + 2 ( c 2 - ^ 



1 \ 2 2 , «X 
- s m 2 H — 

c-J ~ ///f r 



(B.23) 



x 



\n(m w + m^ + m 2 ) + j (5 - 4c 2 ) m 2 ^ - — ^- (m^c 2 + m 2 s 2 ) | ln(m^ + m 2 ) 
(m| - m 2 ) ln(m z + m 2 ) + m^ |m 2 ^ - c 2 + ^ + m x (\ + -^j j 



m 2 ^ 



+m 



2 i / 5 5 5 



21n | + C2 ) + 



c 2 5 2 s 2 / „ 1 



2 2 c : 



+ ^7 + 



— \ r 1 2 (m 2 ^ - m 2 ) ln(m w + mi) + (m| - m 2 ) ln(m z + mi)| 



9 9 77Zi 



+^4m^ - 2m 2 + t^t^ ln(2m w + mi) + ^ |^2m 

22/1 s2 \ m i /\ 1 A l/ 44l 3 4l \ 

+s m 2 12 lnm 2 + — J + — ^1 + 7^ J + ^y- m 2 lnm 2 + -m x In mi J 



ln(2m^ + mi) 



% 2 2 



{ (^c + gi s) 4 ln(2m ZL + mi) + (^ 2 s - # x c) 4 ln(2m 7L + m x ) 
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+2g% \n(2m WL + mi) + 2(g 2 c + g l s) 2 (g 2 s - g x c) 2 ln(m ZL + ™ 7 l + mi) 



+4:g%gf (s 2 \n(m WL + m Z L + m 2 ) + c 2 \n(m WL + m lL + m 2 )) } 



V 



8tt 2 /3 2 



0^9 1 10 , Aq 2 n f mt A/ /10 5 

x(ln//- - c, + - + - In 2) - y2 2 f 7 w [-j - 11 In 2 



(B.24) 



V 



9 2 y 



I28ir 2 (3 2 



9mJ - m\ - 3m 2 ) (in /i 2 /3 2 - a + ^ - In 4 j + 48mJ In 2 



(B.25) 



167T 2 /? 2 



m 



(( " ^4 " ^ + | - C ' - j) + ™Z) 



11 ,\ . . 31, _ 2n2 11 51 

^j + ^- 5-4c I ln(2m w + m z ) + — In// /3 - — a - —c 2 



251 
"96" 



(B.26) 



11 5 1 . 1 2/ L 1 / \, 51 , 
c + — -4s 2 n2 + -c 2 c- - + - 2 - — he In - 

12 4c 8c 2 4 v 2 J 8 V cV 4 3 



/ 1 23 1 \ \ 

+ i-^j - — + 5c 2 + -c 4 J lnm ff j - m w m WL (l + c) - 2s 2 m 2 VL ln(2m WL ) 



+ {^r?w ~ 2c2m WL^) ln(2m H/L + m z ) + ]^rn 2 WL 

( 



m 



m w 



I I 2 ., 2 2 v'""' "' :> / 

+s <i | m w - 2m WL - 2m 7L + — 2 



rn 



w 



j ln(m WL + m lL + m w ) 

2 



+m 1 i(mn/t - m w ) + 



m WL m 2 L - m* L (m 2 WL -m 2 L ) ln(m 7L + m WL ) 



m w 
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+c 2 < 



m lL — > mzL 



3AV 



Vn 



32n 2 (3 2 



9/i 2 

P 2 



In - c 2 - 2 ln{m 1 (mi + 2m 2 ) } 



l/ 4 = 



64tt 2 /3 2 



fe 4 /^ 1 1 35 \ n -2.2 v 293 



293 1 13 



18c 2 18c 4 



g\n f /10 5 \ n2 (g\s 2 3g\ 2 



27 Vc 2 c 4 



-71- 14 +9y 



3c 2 4 



+ 4, 2 ) I^ + ^ + a) (l + 



+3g 2 Y \ (in fL 2 (3 2 - Cl + |) - ^{ - 4m^ - 2m| - 24m) In 4 
+ (6m^ + 3m| + m 4 + 3m 4 - 12m)) (in /i 2 /3 2 - Cl + } 



32ti 2 (3 2 



\(m 1 + 3m 2 ){5'2 i^rn WL + Am w + m ZL c 2 + 2m z c' 



+9i (ttizls 2 + 2m z s 2 ) + m lL [gls 2 + g 2 c 2 ) } 

+ |fl , 2fl'i(^i - ^2) {2m z sc + (m ZL - m 7L )sc} + Ag\m WL (m w + m z c 2 ^ 



J rg\m w (4m 7t s 2 + %m w + |m z c 2 + 4m ZL c 2 ) + 3A (m 1 m 2 + |m 2 + fm^ 
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B.3 SU(2)-Higgs model 

In this section, for completeness, the potential of the SU(2)-Higgs model to order g A , A 2 is 
presented. It is readily obtained from the standard model results of section [0T2| by performing 
the limit g%, gy — > and setting the number of families rif to zero. With the masses 



m 2 . 



X(p2 ~ V + 120 2 ( 6A + 5^) m ' = m 2 + 2Av? 2 



ml = I^V + I 1 



™ 2 = 



the final result is given by the sum of the following contributions: 



(B.30) 



1/1, 3 , 



+ * 1 
4 V 12n8 



mi 3 + 3 m 2 3 + 6 m 3 + 3 m L 3 



(B.31) 



zg 2 



1287T 2 /? 2 



2m 2 (in ^ - ln/i 2 /5 2 - ^ci + ^c 2 + ^ j + m 2 (mi + m 2 ) 



Af m l + 3m A (_ 4In £ + ln/i 2 /5 2 -c^ - i| 
2 v 7 \ 3 J m 



(B.32) 



mi — m 2 ) 2 (m 1 + m 2 ) — m(mi + 3m 2 



1 2 

- m? — m 2 . ) ln(mi + m 2 ) + (m 2 — 4m 2 ) ln(2m 2 + 

m 2 v J v ' 



m 



H - |m 4 — 2 (m 2 + m|) m 2 + (m 2 — m 2 ) | In (mi + m 2 + m) 



3g 2 



64ir 2 8 2 



'5 8 5 



o i t» /_/ u , ,\ m>i j- 

m - In — + -c 2 + ln(2mL + mi) + - — - In mi + -m\m 
V4 9/x 8 / 4m 2 2 



(B.33) 



•- — ^-(m 2 — m 2 ) 2 ln(m + mi) + f2m 2 — m 2 + — ^ j ln(2m + mi) + ^m 2 



14 



16tt 2 /3 2 



m 2 / I I 51 145 

— (31 ln/i 2 /? 2 - 66 lnm + 39 In 3 - —c x - —c 2 - — - 102 In 8 



Sm^m + + (^2 m2 — 6m 2 ^j ln(2m^ + m) 



(B.34) 
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3AV 



V„ 



32tt 2 /5 2 



9/i 2 

In—- - c 2 - 21n{mi(mi + 2m 2 )} 



(B.35) 



V, 



m 



(B.36) 



V 2 



32tt 2 /3 2 



■(mi + 3m 2 )(3m£ + 6m) + g 2 m(8m + Ylmi) 



1 5 
+3A(m!m 2 + -m\ + -m^) 



(B.37) 
(B.38) 
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Appendix C 

Formulae for renormalization at T = 



C.l Abelian Higgs model 

The self-energy of the vector field (in Minkowski space) has the structure 

iW{q 2 ) = i (u(q 2 )g^ + IL L q»q v ) . (C.l) 

Here the momentum independent part of n(g 2 ) is just the correction to the mass, and Ui,(q 2 ) 
is irrelevant in Landau gauge because the propagator is transverse. The complete one-loop 
result without one-particle reducible tadpoles reads 

I% 2 ) = l e 2 A(m 2 ) + 



3-2e 



' m 2 — m 2 



l) A{m%) + - l) A(m 2 ) 



(C.2) 



+ 



1 - e )8m z + 



(m 2 - m 2 - q 2 ) 2 



Q 



B(q 2 : ,m 2 ,m 2 ^ 



Similarly, the one-loop self-energy of the Higgs particle is given by 

ie 2 
2m 2 



3zAA(m 2 ) + 3 - 26 - ^ 1 ie 2 A{m 2 ) + ^—{m% - q*)B(q 2 , 0, 0) 



(C.3) 



(3 - 2e)2m 2 - 2q 2 + 



2m? 



B(q 2 , m 2 , m 2 ) + 18iX 2 v 2 B(q 2 , m 2 , mj) 



Here A(q 2 ) and B(q 2 , wig, m\) are the usual one- and two-point functions, defined in appendix 
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C.2 Standard model 



The Landau gauge self-energy of the Higgs particle is given by 



li^q 2 ) = 3i\A(M 2 H ) + UV 18B(q 2 , M 2 H , M 2 H ) + 6B(q 2 , 0, 0) 



(C.4) 



• 4 2 

%g\v 2 



1 A(M 2 W )+ S-2e-^ r + -^-)B( q 2 ,M 2 w ,M 2 w ) 



2M W 



M w 4M W/ 



< w 



w J 



4M 4 , 



c 4 



«02 



-2 + 2e + 



<? 2 \ 



M 2 W ) 



A(M\ 



w) ~ j^W , M w> 0) + J^Btf, 0, 0) ' 



L w 



• —A M w -> M z 



- 3^y [2A(m t 2 ) + (4m 2 - g 2 )5(g 2 , m 2 , m; 



Here Mw , and m t are the zero temperature masses of the vector bosons and the top 
quark and c is the cosine of the Weinberg angle. 

The fermion self energy can be written in the form 



(C.5) 



The 75-contribution is irrelevant for the mass counterterm 5m t = S^(m 2 )+m t Ey(m 2 ), which 
is obtained from 



4K 2 ) 



■^[5(m 2 ,m 2 ,M 2 ) -5(m 2 ,m 2 ,0)] - ^/*(0, m?) 



(C.6) 



»2 



4 S 2 / S (0,m 2 )-^- 3c- S > ' •-■ 



Is(M 2 z ,m 2 t ) 



4K) 



4m? 



2A(m 2 ) - A(M£) + (M 2 - 2m 2 t )B(m 2 t , m 2 t , M H ) - 2m 2 t B(m 2 t , m 2 t , 0) 



-m 2 B(m 2 , 0,0) 



--(? 2 / y (0,m 2 ) 



(C.7) 



^(M 2 ,, 0) + ^ J V (0, m 2 ) + 1 1^ + 1 ( 3C - ^ ) J J^(M 2 , m 2 ) 
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Here the following short forms for typical self-energy contributions have been used: 



Isim^m^ = -imi (3 - 2e)B{m u m ,m 1 ) 



(C.8) 



2m 2 m\ 



{(2-2e)m 2 + m\ — m 2 t } A(m 2 ) - (2 - 2e)m 2 A(m 2 1 ) (C.9) 



+ ■ 



2m 2 m 2 (l - 2e) + (3 - 2e)m 2 (m\ -m\- m t 2 ) 



+ (m\ — ml — m'fj ^jB(m 2 , ml, ml) — (m 2 — mfj B(m 2 ,0,mf) 
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